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Abstract 

We present the theory of nonzero temperature (T) spin dynamics and trans- 
port in one-dimensional Heisenberg antiferromagnets with an energy gap A. 
For T < A, we develop a semiclassical picture of thermally excited parti- 
cles. Multiple inelastic collisions between the particles are crucial, and are 
described by a two-particle 5-matrix which has a super-universal form at low 
momenta. This is established by computations on the 0(3) cr-model, and 
strong and weak coupling expansions (the latter using a Majorana fermion 
representation) for the two-leg S = 1/2 Heisenberg antiferromagnetic ladder. 
As an aside, we note that the strong-coupling calculation reveals a S = 1, two 
particle bound state which leads to the presence of a second peak in the T = 
inelastic neutron scattering (INS) cross-section for a range of values of momen- 
tum transfer. We obtain exact, or numerically exact, universal expressions for 
the thermal broadening of the quasi-particle peak in the INS cross-section, for 
the magnetization transport, and for the field dependence of the NMR relax- 
ation rate \/T\ of the effective semiclassical model: these are expected to be 
asymptotically exact for the quantum antiferromagnets. The results for 1/T\ 
are compared with the experimental findings of Takigawa et.al. and the agree- 
ment is quite good. In the regime A < T < (a typical microscopic exchange) 
we argue that a complementary description in terms of semiclassical waves 
applies, and give some exact results for the thermodynamics and dynamics. 
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I. INTRODUCTION 



For more than a decade now, much effort has been devoted to understanding the prop- 
erties of a variety of insulating one-dimensional Heisenberg antiferromagnets. By now, the 
basic facts about these systems are very well established: Heisenberg antiferromagnetic 
(HAF) chains with integer spins at each site exhibit a gap in their excitation spectrum while 
those with half- integer spins are gapless Among the spin-1/2 ladder compounds, 

those with an even number of legs exhibit a gap just like the integer spin chains |5],f||| while 
ladders with an odd number of legs are gapless analogous to the half-integer chains || . 

Theoretically, the universal low-energy properties of the gapped systems are well de- 
scribed by the one dimensional quantum 0(3) non-linear a-model (NLaM) without any 
topological term 0,03 • A lot is known exactly about this field theory g,|TI],|n|] and this is 
directly useful in understanding the gapped systems. The spectrum of the cr-model consists 
of a triplet of massive spin-1 particles as the lowest energy excitations followed by multipar- 
ticle continua with no bound states. Many zero temperature (T) properties of the gapped 
systems, including low frequency dynamic correlations, can be explained using the exact in- 
formation available on the cr-model ||12||. On the other hand, until very recently exact results 



for T > were restricted to static, thermodynamic properties |T3| while many experimental 
observables (such as the inelastic neutron scattering (INS) crosssection and NMR relaxation 
rates) directly probe dynamical correlations at non-zero temperature. 

(The universal low-energy properties of the gapless systems have been treated via a 
mapping to a certain critical field theory |Hj . In contrast to the NLaM, powerful techniques 
that exploit the conformal invariance of the theory can be used to determine exactly for 
T > some dynamical correlators that are directly probed by NMR experiments |15,16,17| . 



Similar methods have been used to obtain results for T > on static properties as well [18 



Transport properties have also been studied recently p9| , |20|j2T[| , with results that are quite 
different from those we shall obtain here for gapful systems.) 

This paper shall deal exclusively with the T > dynamical properties of gapped Heisen- 
berg spin chains. A portion of our results have appeared earlier in a short report ? where 
we presented them in the context of the continuum NLcxM, but did not fully discuss their 
range of applicability. Here we shall take a more general point of view of working directly 
with lattice Heisenberg antiferromagnets. The main, and essentially only, requirement on 
the spin chain being studied is that it have an energy gap and that its low-lying excitations 
consist of a triplet of spin-1 particles with the dispersion 

e(k) = A+ C ^ + 0(k*)... (1.1) 

Here k is being measured from an antiferromagnetic wavevector Q = n/a (k = q — Q, and a 
is the lattice spacing), and we have introduced a velocity c to parametrize the mass of the 
particles as A/c 2 . This is in keeping with the 'relativistic' spectrum of the 0(3) cr-model 
e(k) = (A 2 + c 2 k 2 y/ 2 , although most of our results will not rely on this relativistic form. 
Gapped spin chains with a spontaneously broken translational symmetry (spin-Peierls order) 
can have spin-1/2 particle excitations: we shall not deal with this case explicitly, although 
we believe most of our results can also be extended to these systems. 
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The energy gap A is an important energy scale which shall play a central role in our 
analysis. The majority of our results will be in the regime T C A (we shall use units with 
h = k B = 1 throughout) which we now discuss. 

In this regime, there is a dilute gas of excited particles present, and their motion 
and collisions dominate the dynamical properties we study. In particular their spacing 
~ ce A/ ' T /(AT) 1/ ' 2 is much larger than their thermal de-Broglie wavelength ~ cj (AT) 1 / 2 . As 
argued in Ref , these particles can be treated classically except when two of them collide. 
Such two-particle collisions need to be treated quantum-mechanically and are described by 
an 5-matrix, which is, in general, a complicated function of the particle momenta and spin 
orientations. Conservation of total momentum and energy implies that momenta before 
and after collisions have to be the same, and 0(3) invariance and unitarity impose further 
constraints, but a fairly complex structure is still permitted-we will see some explicit exam- 
ples in this paper. However, the r.m.s. thermal velocity of a particle vt = (T/ A) l l 2 c — > 
as T/A — > and thus we need the 5-matrix only in the limit of vanishing incoming (and 
outgoing) momenta. One of the central ingredients in our computations will be our claim 
that in this limit, all of the complexity disappears, and the 5-matrix has a super-universal 
form; for the scattering event shown in Fig |l] we have (here (ii — x,y,z are the three possible 
values of the 0(3) spin label) 

h; k[, k' 2 ) = (-1)6^6^(2^6^ - k[)5(k 2 - k' 2 ) ; (1.2) 

Notice especially the opposite pairing of momentum and spin labels: crudely speaking, the 
momenta go "through" the collision, while the spins "bounce off" — this dichotomy will be 
crucial to our considerations. We dub this limiting value of the iS-matrix 'super-universal' as 
it requires only that the lowest lying excitations above the gap satisfy ( |1 . 1|) at low momenta. 
The value, however, does not depend on parameters such as c and A. Moreover, we expect 
this limiting result to hold even at the lattice level for generic microscopic models of one- 
dimensional antiferromagnets with massive spin one excitations; we shall see one explicit 
example that bears out this expectation later in the paper. 

With this simple form of the iS-matrix in hand, we will use the semiclassical techniques 
of Ref ]23[ to analyze dynamical properties of spin fluctuations near q = and q = Q in 
terms of the motion of the dilute gas of quasiparticles. 

We will begin by discussing the properties of the iS-matrix for two-particle scattering in 
the limit of low momenta in Section f|. 111 Section |II A| we consider the S matrix for the 
0(3) non-linear sigma model. This has been computed for all momenta by Zamolodchikov 
and Zamolodchikov, and we shall show that the zero momentum limit does indeed satisfy 

However, the cr-model is a continuum theory; it would be much more satisfactory to 
be able to directly see that (|1.2| ) holds for some specific microscopic model, and explicitly 
verify that lattice effects do not affect the simple structure of this limit. One such model 
is the S = 1/2, two-leg Heisenberg antiferromagnetic ladder with inter-chain exchange Jii 
and intra-chain exchange J±. The properties of the ladder can be analysed using a 'strong- 
coupling' expansion |24|] in powers of J\\ / J± for the microscopic lattice Hamiltonian of the 
system. In Section [11 B| we shall explicitly verify ( |1.2| ) for vanishing velocities in this lattice 
model within this strong-coupling expansion. Parenthetically, we note that our strong- 
coupling analysis also allows us to make predictions about interesting features in the T = 
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dynamic structure factor S(q, uj) which are specific to the system considered. In particular, 
we find that, to second order in J\\/J±, a two-particle S = 1 bound state gives rise to a 
second peak (in addition to the usual peak coming from the stable single particle excitations 
of the system) in S(q, uj) for a range of values of q around Q. This should be of relevence 
to inelastic neutron scattering (INS) experiments on the ladder compounds and it is hoped 
that they experimentally verify the existence of this effect. 

In Section [11 (J| we study the complementary 'weak-coupling' expansion in powers of 
J±/J\\ for the two- leg ladder. As was shown in Ref pB |, this expansion leads to a description 
of the low-energy, long-distance properties of the ladder in terms of an effective field theory 
of a triplet of massive Majorana fermions. The Hamiltonian for the Majorana fermions also 
has a four-fermion coupling which has generally been ignored in previous treatments. In 
the absence of this scattering, the Majorana fermions are free, and the resulting iS-matrix 
does not obey (|1.2| ). In this paper, we consider the effect of the four-fermion coupling 
in perturbation theory. We show that this expansion suffers from severe infra-red problems 
which have to be resolved by an infinite-order resummation. The structure of the divergences 
is very similar to those also present in the large- N expansion of the cr-model above, and we 
find that the resulting resummed 5-matrix of the Majorana fermions does indeed obey the 
analogue of ( |1.2| ). So neglecting the four-fermion coupling | ]26[| , (or even treating it in an 
unresummed manner at finite order in perturbation theory) is a very bad approximation at 
low momenta, and we expect that corresponding divergences in the perturbative evaluation 
of the spin-spin correlation function invalidate the dynamical results of Ref |26j at low T. 

In Section |T| we shall turn to a discussion of the dynamical properties in the regime 
0<T<A. Our results apply universally to all gapped one-dimensional antiferromagnetic 
systems with spin one quasiparticles; indeed they rely only on the dispersion (|1.1| ) and the 
iS-matrix in ( p..2[ ). All our results will be expressed solely in terms of the parameters c and 
A, the temperature, T and the external field H. 

In Section |111 A| we study the dynamics of the staggered component (with wavevector q 
close to Q) of the fluctuations in the spin density. More precisely, we study the dynamical 
structure factor S(q,uj) for q close to Q. Apart from some overall factors, this directly gives 
the INS crosssection at the corresponding values of momentum and energy transfer. At 
T = 0, the dynamical structure factor has a sharply defined 5-function peak at uj = e(q — Q) 
for q near Q. This peak can be thought of as arising from the ballistic propogation of the 
stable quasiparticle of the system. At non-zero temperatures, the peak broadens as the 
quasiparticle suffers collisions with other thermally excited particles. The main objective of 
Section [III A| is to describe the precise lineshape of the quasiparticle peak in the dynamic 
structure factor for T > 0. 

In the cx-model approach, the staggered components of the spin density are represented 
by the antiferromagnetic order parameter field n. We will use the semiclassical method of 
Ref- |23[ to calculate the space and time dependent 2-point correlation function of the n 
field for T > 0. This allows us to calculate the thermal broadening of the single particle 
peak in the dynamical structure factor S(q,uj) for wavevectors q near Q. In particular, we 
find that the dynamic structure factor in the immediate vicinity of the quasiparticle peak 
at (q = Q, uj = A) may be written in a reduced scaling form as 
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S(q,u) 



AcL t I oj — e(k) 
T^aT 



LI 1 



;i.3) 



where k = q — Q, A is the (non-universal) quasiparticle amplitude of the spin one excitations 
of the system, L t = A/vre A//T /3T is the typical time spent by a thermally excited quasiparticle 
between collisions with other particles, and <£> is a completely universal function that we 
determine numerically in this paper. Notice that temperature enters this scaling form only 
through L t . We claim, though this is not rigorously established, that these results for 
the broadening are asymptotically exact for T C A: all corrections to the line- width are 
expected to be suppressed by positive powers of T/A. Some evidence for the exactness of 
our results emerges from consideration of simpler systems where exact results for the line- 
broadening are available from the quantum inverse scattering method |j27 |; as we shall see in 
Section [III A| , our semiclassical results are in perfect agreement [23[] with these. It is hoped 
that experimental studies of the temperature dependence of the INS cross-section in this 
regime will confirm these results, particularly the simple scaling form (ITT 



In Section [III B| we turn to the correlations of the conserved magnetization density, 
or dynamic fluctuations near q = 0, for T < A. Unlike the staggered case, the overall 
magnitude of the magnetization density fluctuations is universal and given by Tx u , where 
Xu is the uniform susceptibility of the system (the non-universal overall scale of the staggered 
component is reflected for instance by the presence of the overall constant A in ( |1.3| )). In 
this temperature regime, we have the well-known result for Xu [|13||: 



Xu 
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We shall study the dynamics of the magnetization density in Section [111B| [|22|. We shall 



show that the long-time correlations of an effective semiclassical model are characterized by 
spin diffusion, and obtain the following result for its low T spin-diffusion constant D s : 



3A 



;i.s) 



Using the Einstein relation for the spin conductivity a s = D s x U i we obtain from (|1.4|) and 



cf 2 \V a 
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Notice that the exponentially large factor e A//T has dropped out, and a s diverges with an 
inverse square-root power in T as T — > 0. The semiclassical model is possesses an infinite 
number of local conservation laws: in Appendix A, we discuss how the existence of spin 
diffusion can be compatible with these local conservation laws. However, these results do 
not rigorously allow us to conclude that the ultimate long-time correlations of the underlying 
gapped quantum spin chain are diffusive. This has to do with a subtle question of order 
of limits: we computed the S matrix ( |1.2j ) in the limit T/A — > 0, and then used it to 
evaluate the long-time limit of correlations of the magnetization, whereas in reality the 
limits should be taken in the opposite order. What we can claim is that our results will 
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apply for all times upto a time scale which is larger than the collision time L t by a factor 
which diverges with a positive power of A/T as T — > 0; there is a substantial time window 
in this regime where we have established that the spin correlations are diffusive. For the 
generic gapped quantum spin chain, we can reasonably expect that the ultimate long-time 
correlations are indeed diffusive, and the only consequences of the omitted terms in the S 
matrix are subdominant corrections to the value of D s in ( |1.5|) which are suppressed by 
powers of T/A. For the continuum NLoM with a relativistically invariant regularization 
(this is unphysical for any experimental application), the issue is a little more subtle: this 
model does possess additional non-local conserved quantities |28| , but we consider it unlikely 
that these will modify the long time limit |^9| . On the experimental side, however, diffusive 
behaviour of the magnetization density has already been convincingly demonstrated in the 
S = 1 one-dimensional antiferromagnet AgVP2S6 by the NMR experiments of Takigawa 
et.al. 0. 



As has been argued earlier |pi|, the dynamic fluctuations near q = provide the dominant 



contribution to the NMR relaxation rate 1/Ti for T < A. Thus, knowing the space and 
time dependent two-point correlation function of the conserved magnetization density, we 
are able to compute the field and temperature dependent 1/T\ in this regime. We shall see 
that the overall scale of 1/Tj is set by the ratio Tx u / VDs- As was pointed out to us by M. 



Takigawa |22|, this immediately leads to an activation gap for 1/Ti given as 



This difference between the activation gaps for \u and 1/Ti appears to clear-up puzzling 
discrepencies in the experimental literature |30,32,5| for the value of the energy gap in these 



systems obtained from Knight-shift susceptibilty measurements on the one hand, and 1/Ti 
NMR relaxation rates on the other; a systematic tabulation of the activation gaps for a large 
number of gapped spin chains [33| does indeed show a trend consistent with ( |1 . 7| ) . The crucial 



factor of 3/2 clearly arises from the exponential divergence in D s . This diffusive behaviour 
we find arises entirely from intrinsic inelastic scattering between the quasi-particles. In real 
systems there will also be contributions from elastic scattering off inhomogenities which will 
eventually saturate the divergence of D s as T — > 0. However, because of the strong spin 
scattering implied by ( |1.2| ) the effects of inelastic scattering is particularly strong in d = 1, 
and can easily dominate inhomogenities in clean samples. 

We will give a detailed account of the calculations leading up to our expression for 1/Ti 
(some details on the method used are relegated to Appendix C) and then go on to compare 
the theoretical predictions for the field dependence of 1/Ti with the extensive experimental 
data of Takigawa et.al. |3(J. We will see that our results {without any adjustable parameters, 
except for a field independent background rate) agree with the data extremely well for a 
range of intermediate temperatures. At the lowest temperatures for which data is available, 
the quality of the fit deteriorates significantly and the 1 / yH divergence predicted at small 
fields seems to get cutoff, presumably by some spin-dissipation mechanism present in the real 
system. At the present time, we are unable to incorporate this dissipation in any serious 
way in our approach. However, following Ref JJIJ, we can phenomenologically introduce 
some spin-dissipation in our results for the long-time limit of the autocorrelation function 
and obtain a corresponding expression for the field dependence of l/Ti. This allows us to fit 
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the data at the lowest temperatures with a phenomenological form that has one additional 
adjustible parameter corresponding to the spin-dissipation rate. We also present results for 
the temperature dependence of this effective rate. 

Finally, in Section [IV| we will turn to the regime T ^> A. We will do this in the 
context of the continuum 0(3) a-model only. Any continuum theory is applicable to real 
lattice experiments only below some energy scale, and a natural choice for this energy scale 
is a typical exchange constant J. So more specifically, we shall be studying the regime 
A C T C J. For T > J we expect the spins to behave independently, and the system 
exhibits a Curie susceptibility. It is an open question whether the window of temperatures 
A < T <C J with universal behavior exists at all in any given system, and the answer 
will surely depend upon details of the microscopies. It is unlikely to be present for 5=1 
spin chains, but appears quite possible for S = 2 spin chains [j3H}| . The static properties 



of this regime were first studied by Joliceour and Golinelli |35[ using the N = oo limit 
of the 0(N) a model of Ref ||36| . We shall present here an exact treatment of static and 
dynamic properties for the case of general 0(N); the numerical values of the N = 3 static 
results are significantly different from the earlier N = oo results. We shall show that the 
antiferromagnetic correlations decay with a correlation length £, which to leading logarithms 
in A/T is given at iV = 3 exactly by 

s 2nT V A / ' K ' 

where 7 is Euler's constant. We also obtain the exact uniform susceptibility 

(notice the argument of the logarithm differs slightly from (|l.<j|)). It is interesting to compare 



the two asymptotic results (|1.4j) and (|1.9|) , and we have done that in Fig |]. It is reassuring 
to find that the two results are quite compatible for T m A. This suggests that one of 
either the T <C A or T > A asymptotics are always appropriate. We shall also consider 
the nature of spin transport in the A C T C J regime, and show that it is related to 
transport in a certain classical statistical problem of deterministic non-linear waves. We 
have not established whether spin diffusion exists or not in this classical problem; if the 
correlations were diffusive, however, we are able to precisely predict the T dependence of 
the spin diffusivity: 

D = B- — n (1 10) 



Here B is an undetermined universal number, and £(T), Xu(T) are given in ( |1.8| . |1.9| ). 

Notice the complementarity in the two T regimes discussed above: the description for 
T C A was in terms of semiclassical particles, while that for A <C T <C J is in terms of 
semiclassical non-linear waves. 
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II. ZERO TEMPERATURE PROPERTIES 



The primary purpose of this section will be to establish the iS-matrix by a variety of 
methods. We will begin in Section [II A| by using the relativistic 0(3) a model. In Section [HB 
we will consider the strong-coupling expansion of the two-leg ladder in powers of J«/ J±. This 
section will also present supplementary results on some interesting features in the T = INS 
cross-section of the strongly coupled ladder arising from the presence of a S = 1, two-particle 
bound state in its spectrum. Finally, Section [11 C| will consider the complementary J±/J\\ 
expansion. 



A. (9(3) cr model 

Let us begin with a brief review of the cr-model as an effective field theory for the 
low-energy properties of the gapped systems (for a more extensive discussion see |37| and 
references therein). The imaginary time (r) action of the cr-model is 



c r x l T \ 1 
£ = — / drdx (d x n a ) 2 + —(d T n a - ie^Han 
2g Jo L cr 



) 2 



(2.1) 



where x is the spatial co-ordinate, a,/3, 7 = 1,2,3 are 0(3) vector indices over which there 
is an implied summation, e a ^ is the totally antisymmetric tensor, c is a velocity, H a is 
an external magnetic field (we have absorbed a factor of the electronic magnetic moment, 
g e fiB, into the definition of the field H.), and the partition function is obtained by integrating 
over the unit vector field n a (x, r), with n^(x, r) = 1. The dimensionless coupling constant 
g is determined by the underlying lattice antiferromagnet at the momentum scale A ~ 
inverse lattice spacing to be g ~ 1/S where S is the spin at each site in the original lattice 
system. The cr-model is used to make statements about physics at length scales 3> A -1 
and time scales 3> (cA) _1 ; this physics is universally characterized by the dimensionful 
parameters c, H, T, and A, the energy gap at T = H = 0. Though the magnitude of A 
is determined by non-universal lattice scale quantities (A ~ cAe^ 2n ^ 9 for small g), the long 
distance physics of the cr-model depends on these lattice scale effects only through the value 
of A, and has no direct dependence on g or A. Also, the energy-momentum dispersion of 
the stable particle-like excitation of this model is given by e(k) = VA 2 + k 2 c 2 , and there is 
a triplet of them. The conserved density of this model corresponding to its 0(3) symmetry 
is the magnetization density M a (x,r) = 5L/5H a (x,r). In the Hamiltonian formalism, this 
is represented by the operator M a (x). 

Finally, to make contact with the lattice antiferromagnet, we must have a prescription for 
representing the spin-density operator s a (x) of the lattice system in terms of the operators 
of the cr-model. It is most convenient to do this in terms of Fourier components. We have: 

s a (k + Q) oc n a (k) 

(recall that Q = ir/a) and 

s a {q) = M a (q) 
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for \q\, \k\ much smaller than some microscopically determined scale ~ A. The missing 
proportionality constant in the first relation is non-universal and related to the magnitude 
of the spin at each site in the original lattice system. Thus, the a-model allows us to represent 
spin fluctuations near q = Q (these being the low-energy degrees of freedom) and near q = 0. 
This is of course because the q = component of the spin density is the conserved charge 
corresponding to the 0(3) symmetry of the system, and as such must be included in any 
description of the slow modes. 

The exact S matrix of the collision of two particles in the cr-model was computed in 
a seminal paper of Zamolodchikov and Zamolodchikov. For the scattering event shown in 
Fig [L]it is (recall /ij — x,y,z are the three possible values of the 0(3) spin label): 

k 2 ; k[, k' 2 ) = (2tt) 2 5(A; 1 - k[)5{h - k' 2 ) M^ 1M2 V 



+^2(#)W 1 ( W + ^WW^'Wi ( 2 - 2 ) 



where 6 = 6\ — 6 2 is the 'rapidity', ki = (A/c)sinh#j for i = 1,2, and 0(3) invariance 
guarantees a total lack of H dependence in the result. The functions a in (|2.2| ) are 

<T 2 (0) 



2md 


[6 + m){6 - 


2m) 


6(6 - TTI 


) 


(6 + m)(6 - 


2vri) 


-2ni(6 - 


ni) 


(6 + TTi)(6 - 


2iti) 



= T7T7 -In o ( 2 - 3 ) 
(6 + m)(6 — 2m) 

Now notice the structure of the limit 6 — > which is important for our purposes in the 
region T <C A: we find that 012(0 — >• 0) = 0, while a^(6 — > 0) = —1. This establishes the 
key result ( |1.2j ) for this continuum model. 



B. Strongly-coupled two-leg ladders 

In this section we concentrate on the properties of a particular model system, the spin- 
1/2, 2- leg Heisenberg antiferromagnetic ladder, to which the low-energy phenomenology of 
the preceeding section is expected to apply. The Hamiltonian of the system may be written 

as 

n = £ s T (i) ■ s n {i) + • s T {i + 1) + s n {i) ■ ^ + 1)) . (2.4) 

i i 

Here, the Sj(i) and Sjj(i) are spin-1/2 operators at site i along the two chains that make up 
the ladder, g is a dimensionless coupling constant equal to the ratio of the antiferromagnetic 
bond along the individual chains Jy to the bond along the rungs of the ladder Jj_ and we have 
set the bond strength J± along the rungs to be unity; this defines our unit of energy. We will 
analyze this model in the limit of small g; this 'strong-coupling' expansion EM is expected 



to be qualitatively correct for all g. For g = 0, we just have a system of isolated rungs with 
the two spins on each rung coupled antiferromagnetically. The ground state is a product 
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state with each rung in a singlet state. The lowest lying excited states form a degenerate 
manifold with precisely one rung promoted to the triplet state. Perturbative corrections in g 
would presumably make this triplet 'particle' hop around producing a single-particle band of 
triplet excitations as the lowest lying excited states. Thus we expect that our perturbative 
analysis will be most conveniently performed in a representation that directly describes the 
state of individual rungs of the ladder. With this in mind, we switch to the 'bond-operator' 
formalism introduced in Ref [S8|. Following Ref [ fj8|| , we write the spin operators as: 



S?(i) = \ (s\i)t a {i) + ft (i>(i) - ie a ^(i)t 7 (i)) , (2.5) 
Sfi® = ~ (-sHi)t a (i) - £(<)*(<) - ie^(«)*7(0) , (2-6) 

where a, /3, and 7 are vector indices taking the values x,y,z, repeated indices are summed 
over, and e is the totally antisymmetric tensor. s'(i) and v a (i) are respectively creation 
operators for singlet and triplet bosons at site % (in the previous section we had used s(x) to 
denote the spin density of the lattice system; here we shall use a(x) to denote the same and 
reserve s for the singlet boson operator). The restriction that physical states on a rung are 
either singlets or triplets leads to the following constraint on the boson occupation numbers 
at each site: 



'(*)*(*) +£(*%(*) = 1 



The spin density is given by 



a a (i) = -160^^(1)^(1) . 

It is also convenient to define 

4> a {i) = s f («%(*) + *«(*)*(»") • 
The Hamiltonian in terms of these operators is given as 

H = H + V ; (2.7) 



where 



Ho = E (-^(iMO + , (2.8) 



and 



V = I E (0a(«>a(i + 1) - <Ta(i)<Ta(i + 1)) • (2-9) 

In this representation, the ground state for g = is just the state with every site occupied 
by a singlet boson. To zeroth order in g, the lowest excited states form a degenerate manifold 
with a triplet boson (of arbitrary polarization) replacing the singlet particle at precisely one 
site. Higher excited states also form degenerate manifolds labelled by the number of singlet 
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particles that are replaced by triplet bosons. In what follows, we will describe states by the 
number (which can only be zero or one) and polarization of the triplet particles at each site, 
the singlet occupation numbers being determined by the constraint. Thus we will loosely 
refer to the state with no triplet particles as the 'vacuum'. At this order in g, the physical 
particle-like excitation of the system is created at site i by the action of 4> a {i) on the vacuum, 
and thus coincides with the bare triplet particle. In general at higher orders in g, we expect 
that the physical single particle states of the system will contain an admixture of states 
with more than one bare particle present. Similarly, the physical vacuum will also have a 
component with non-zero bare particle number. 

In fact, it is quite convenient to make a canonical transformation (determined order by 
order in g) to an auxiliary problem in which the physical particle states do not contain 
admixtures of states with different bare particle number. The Hamiltonian of the auxiliary 
problem is related to the original one by a similarity transformation. The energy eigenvalues 
obtained in this manner of course give the energy levels of the original Hamiltonian. However, 
to recover the corresponding wavef unctions, one has to undo the effects of the canonical 
transformation. We will use this convenient formulation of perturbation theory below as we 
discuss the strong-coupling expansion. 

The auxiliary Hamiltonian in this approach is given by 

U = e lW He~ lW , (2.10) 

where W is the hermitian operator that generates our canonical transformation. We choose 
W to meet the following criteria: 

• The matrix elements of 7i between states with differing numbers of bare particles 
should be zero to a given order in g. Note that this implies that the elementary excita- 
tions of the auxiliary Hamiltonian are just the bare particles. However their dynamics, 
and their mutual interactions (in multiparticle sectors) are determined by the restric- 
tion of H to the appropriate subspace of definite particle number. This restriction gives 
the corresponding energy levels of the original Hamiltonian 7i correct to that order in 
g. This then serves as our effective Hamiltonian for the corresponding sector of the 
original problem. 

• This does not completely specify W . We therefore also require that W have zero 
matrix elements to a given order in g between any two states populated by the same 
number of bare particles. 

These criteria fix W uniquely order by order in g and in general we have an expansion for 
W that reads: W = g (Wi + gW 2 + ■ ■ ■ ) . It is quite straightforward to use this procedure 
to generate an expansion in g for the effective Hamiltonians in the one and two particle sector 
of the original problem (the 'effective Hamiltonian' in the physical vacuum sector is just a 
constant equal to the ground state energy calculated to the relevant order in g). Solving 
for the eigenstates and eigenvalues of these effective Hamiltonians is just a simple exercise 
in elementary quantum mechanics. If the eigenstates of the spin-ladder are of interest (as 
they will be when we calculate S(q, uj) perturbatively) , we will have to obtain them from 
the eigenstates \ip) of the effective Hamiltonian using 
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m ph y Si cal = e- m \ij) . (2.11) 

After this preamble, we turn to the actual calculations. As we have mentioned earlier, the 
scattering matrix in the low-energy limit is a crucial input to the semiclassical calculations at 
non-zero temperature, and it is therefore interesting to have results for it in our microscopic 
model. So, to begin with, let us look at the two-particle sector and work out the scattering 
properties of the physical particles. 

First we need to find the effective Hamiltonian for the two-particle sector. To first order 
in g, this is just given by the restriction of Ho + V to the two particle subspace. Instead of 
introducing a lot of cumbersome notation to write this down, we will just list the amplitudes 
of the various processes that are allowed in this two body problem: 

• Each particle can hop one site to the left or the right with amplitude g/2 except when 
the neighbouring site in question is occupied by the other particle. 

• When the two particles are at neighbouring sites, there is a non-zero amplitude for 
spin rotation. Consider the state \i,ai, (i + l),^) which has one particle at i with 
polarization a\ (which can be any one of x,y,z) and another particle at i + 1 with polar- 
ization a 2 . The amplitude to make a transition from this to the state \i, j3\] (i + 1), /3 2 ) 

i s ( fi ,e 7/3iai e 7/32«2) • 

To solve for the scattering states of this two-body problem, it is more convenient to work 
in a basis in which we label the spin part of the two particle states by the total angular 
momentum J and the value of its z component J z . The spin rotation amplitude now becomes 
just a J dependent nearest neighbour potential which takes the values g/2, —g/2, and —g 
for J = 2, 1, and respectively. Note that the potential energy is independent of J z as 
one would expect from rotational invariance. It is now quite simple to find the scattering 
eigenstates in each channel. The spatial wavefunction in channel J may be written as 

^( Xl ,x 2 ) = pj{(e lk ^ +ik ^ +rj {k 1 ,k 2 )e ik ^ +lk ^) 6(x 2 - Xl )} , (2.12) 

where Pj is the symmetrizing operator for J = 2, and the antisymmetrizing operator for 
J = 1, and r\j is the reflection coefficient that completely specifies the scattering properties 
of the particles. For the r j we have: 

_e~ ika -2 cos (k cm a/ '2) 
T ' 2 ~ e ika - 2 cos (k cm a/2) ' 

e -ika + 2cos{k cm a/2) 
Tl ~ e ika + 2cos{k cm a/2) ' 

e -ika + cog (k cm a/2) 

T ° ~ e ika + cos (k cm a/ '2) ' 

where k — (ki — k 2 )/2 , k cm — k\ + k 2 , and a is the lattice spacing along the length of either 
of the two chains that make up the ladder system. Note that k\ and k 2 both range over the 
interval (0, 2-k /a). The energy of the scattering state labelled by {k±, k 2 } (the energy of the 
ground state being set to zero) is given by 
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(2.13) 
(2.14) 
(2.15) 



E(ki, k 2 ) = 2 + g cos(fcia) + g cos(/c2a) • 

This is consistent with the first order result |24| for the single particle dispersion rela- 
tion: E(k) = 1 + g cos(ka). 

The next step is to use these results for the reflection coefficients to obtain the 5-matrix 
for this two-body problem in the limit of low velocities. Low velocities imply values of k\ 
and k 2 in the vicinity of the band minimum at n/a, i.e k close to zero and k cm close to 
2ir/a. Both r 2 and r\ have the limiting value — 1 as k — > 0, k cm — > 2ir/a. However, r$ is 
singular in the vicinity of k — 0, k cm = 2-nja\ its value depends on the order in which the 
two limits k — > and k cm -^2-nja are taken. This is somewhat disconcerting as we expect 
a well-defined low- velocity limit which agrees with the predictions of the 0(3) NLcrM field 
theory. 

To identify the source of our problem, let us look more closely at the expression for r . We 
notice that r , considered a function of the complex variable k, has a pole in the upper half- 
plane for a range of values of k cm . This indicates the presence of a bound state in the J = 
channel for the corresponding values of k cm . This bound state hits threshold, i.e its binding 
energy goes to zero, as k cm — ► 2n/a. It is the presence of a bound state at threshold that 
causes the singular behaviour of the reflection coefficient in the limit k — > 0, k cm — > 2n/a. 
Clearly, if there were a range of k cm around 2tt /a for which there was no singlet bound state, 
then we would not have this difficulty. It turns out (as we shall briefly outline later) that 
extending our calculation to the next order in g leads us to precisely this conclusion and 
gives a well-defined limiting value of —1 for r as k — > 0, k cm — > 2% /a. 

This result can now be used to obtain the 5-matrix of our auxiliary two-body problem. 
We are interested, however, in the 5-matrix that describes the scattering of the physical 
particle-like excitations of the spin-ladder. Thankfully, it is quite easy to see that though 
the wavefunctions of the two problems are related by a canonical transformation, the purely 
'off- diagonal' form of W implies that the two are the same at least to first order in g. 
Transforming to the basis used in (|1.2p , we see that the 5-matrix in the low velocity limit is 
indeed given by ( |1 . 2| ) Thus, this super- universal form of the iS-matrix holds for our lattice 
model and lends support to the idea that it is a generally valid consequence of just the 
slow motion of the particles and is in no way dependent on the special properties of the 
continuum cr-model. 

To wind up this part of our discussion, let us now summarize the calculation of the 
reflection coefficients to first order in g. We need to find the effective Hamiltonian of our 
auxiliary two-body problem to second order in g. This involves first working out W\ and 
then using this to obtain the effective two-body Hamiltonian. To 0(g 2 ), we generate in this 
manner an additional next-nearest neighbour hopping term and some additional J dependent 
nearest-neighbour interactions. We skip the details as they are somewhat tedious and not 
particularly illuminating. The J = reflection coefficient (correct to O(g)) obtained in this 
manner is given as: 

e -ika + cog (k cma /2) -gcos {k cm a) (3e~ ika + e ika ) /4 

r ° ~ e ika + cos (k cm a/2) - gcos (k cm a) (3e ika + e~ ika ) ~ji ' ^ 2 ' 1 ^ 

From this, it is easy to see that there is no pole in the upper half fc-plane as long as 
\k cm a — 2n\ < y/Eg. This means that there is no singlet bound state possible in this range of 
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k cm . This is consistent with our expectation that at the very lowest energies, the two-particle 
spectrum should be free of bound states in order to match the predictions of the cr-model. 
Moreover, ( |2.16| ) has a well-defined low- velocity limit of —1 as claimed earlier. 



The foregoing analysis has shown that the two particle sector has a spin S = bound 
state which leads to some interesting threshold singularities for the scattering matrix. Ex- 
amining our expressions for the reflection coefficients, we notice that there is in fact a bound 
state in the S = 1 channel as well (actually, there is also a S = 2 anti-bound state; we will 
not delve further into that aspect of the spectrum here). Now, a S — 1 excited state can 
have observable consequences for the INS cross-section of a system and we might expect to 
see some interesting features in the same as a result of this. 

With this motivation, let us turn to the perturbative calculation of the dynamic structure 
factor at T = 0. We pick a coordinate system in which the two chains that make up our 
ladder are parallel to the x axis and have y coordinates of +d/2 and — d/2 respectively, 
where d is the distance between them (for simplicity, we are assuming here that the rungs 
of the ladder are perpendicular to its legs). The spins along a chain are located at x values 
equal to integer multiples of a. We denote the position of each spin in the x-y plane by R. 
We define P = (q x , q y ) . The T = dynamic structure factor may be written as 



S(P,u) = ^j dt ^H ($o|^(t)%(0)|$o)e-^ (i? - i? ' )+ ^ ; (2.17) 

R R' 



where |$ ) is the ground state of the system, L is the length of each chain, and denotes 
the spin operator at R in the Heisenberg representation. Our strategy is to write down 
the usual spectral representation for ( f2.17| ) and then evaluate it perturbatively. Actually, 
a complete calculation of the second order contribution would involve the eigenstates with 
more than two-particles present; below we will ignore this complication and confine ourselves 
to calculating the contribution of the one and two-particle sectors, correct to the appropriate 
order in g. 

The spin operator at any site is a sum of two terms: a single particle piece coming 
from the (f) a , and a two particle part coming from the spin density operator a a . From the 
structure of the strong-coupling expansion, it is clear that the single particle part does not 
have matrix elements between the ground state and any state in the two particle sector; 
similarly the two-particle piece does not have matrix elements between the ground state and 
any state in the single-particle sector. Thus, keeping only the contributions from the one 
and two particle sectors, we can write to second order in g: 

s(P, = \ E s{<»- Ei) I (®i\M-<i*)\*o)\ 2 sin 2 ( q -f) 

1— particle states \ / 

+ \ E ^-^,)K^k.(-a.)l*o>| a oos a ^) ; (2.18) 

2— particle states \ / 

where |$i) and |$ 2 ) denote one and two particle states respectively, and (f> z {Qx) an d cr z (q x ) 
denote the discrete Fourier transforms of <p z (x) and a z (x) . Let us digress for a moment and 
think in terms of the inelastic neutron scattering cross-section for a process with momentum 
transfer P and energy tranfer a>; this coincides with the dynamic structure factor apart from 
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some geometrical factors. This scattering can of course produce a single spin-one particle 
in the spin system. But there is also a non-zero amplitude for producing a pair of these 
particles close to each other (as is clear from the actual calculations described later). This 



is the origin of the second term in ( pTT 



Now, these two pieces contribute to the structure factor over very distinct intervals along 
the frequency axis. While it is in principle possible to calculate both terms correct to 0(g 2 ), 
we will confine ourselves below to calculating the leading perturbative correction for each 
value of uj. Thus, we will calculate the single particle piece only to first order in g, while 
doing a full second order calculation for the two-particle piece. Below, we give a brief outline 
of the calculation and then discuss our final results. 

To calculate the single particle piece, we first need to determine the ground state and the 
physical one-particle state wavefunctions correct to O(g). This involves using W correct to 
first order to obtain the physical wavefunctions from the wavefunctions of the corresponding 
auxiliary problem (for the one-particle sector, these are just plane waves to all orders in g; 
this follows from translational invariance). A simple calculation then gives the one-particle 
piece as 

S 1 (P, u) = ^ (1 -g cos (q x a)) sin 2 (^y(u-E(q x )) , (2.19) 

where E (q x ) — 1 + g cos (q x a). 

Turning to the two-particle piece, we see that one can actually ignore the distinction 
between the physical 2-particle wavefunction and the wavefunction of the auxiliary two-body 
problem. Moreover, it suffices to consider the auxiliary problem to first order in g. Also, 
since the ground state has spin zero and we are looking at the matrix elements of a vector 
operator, we need to consider only the triplet (J = 1) channel of the auxiliary problem. 
The only subtlety lies in the fact that we need to consider the bound state contribution 
as well as the usual contribution of the scattering states. From ( |2.14| ), we see that this 



bound state exists for n < k cm a < n + ir/3 and for 3n — n/3 < k cm a < 3ir (remember 
k cm ranges from to 4ir/a). Thinking in terms of an inelastic neutron scattering event 
with momentum transfer q x in the fundamental domain (0, 2n/a), we see that this bound 
state can be excited for 7r < q x a < it + it/ 3 and for ir — tt/3 < q x a < rr. In the latter case, 
momentum is conserved modulo a reciprocal lattice vector of 2ir/a. Of course, in addition to 
the bound state contribution there is a background term coming from the scattering states 
in this channel. Again, the two particles can be created in the scattering state either with 
total momentum k cm equal to the momentum transfer q x , or with the two differing by a 
reciprocal lattice vector of 2n/a. 

The actual calculations are quite elementary and we proceed directly to the results for 
the two-particle contributions. The bound state contribution for 2nj3 < q x a < Ati/3 may 
be written as: 



S B {P,u) = ^cos 2 fafj sin 2 (M) (i -4 cos 2 (^)) S (u, - E B (q x )) , (2.20) 



where E B (q x ) = 2 — g(l + 4 cos 2 (q x a/ 2)) / '2 . On the other hand, the scattering states give 
rise to the following background contribution for \uj — 2\ < +2g\ cos(q x a/2)\ : 
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o ( B n 9 2(%d\ ■ 2 fQxa\ y4gcos^g/2) - (a; - 2) 2 

d sc (P, to) = — cos ^— sm — ; — — . 2.21 

V ' ; 7T V 2 / \ 2 J g + 2(u -2) +4gcos 2 (q x a/2) K ' 

Note that for q x a = 2tt/3 or 47r/3, there is a square-root divergence at the lower threshold 
to the continuum in u; these are precisely the values of q x for which the binding energy of 
the triplet bound-state goes to zero. This enhanced scattering can thus be thought of as 
arising from the presence of the triplet bound state at threshold. The salient features of 
these results are summarized in Fig |3] and Fig |]. Fig |3| is a plot of the positions along the 
uj axis of the single particle peak, the bound state peak, and the bottom of the two-particle 
continuum as a function of q x . In Fig |], we show the spectral weight in the single particle 
and bound state peaks as a function of q x . 

Thus, we see that that the existence of a triplet bound state of two elementary spin-one 
excitations leads to some interesting features in the dynamic structure factor. Actually, 
qualitatively similar features, again arising from a triplet bound-state, had been predicted 
earlier in the alternating one-dimensional Heisenberg antiferromangnetic chain. Recent 



INS experiments [0,0 on (VO) 2 P20 7 do indeed see a second sharply defined peak in the 
dynamical structure factor for a range of values of q x . While this compound had been 
previously thought to be a good example of a spin-ladder [|^], more recent work JEJ has 



favoured the alternating chain model |Q and the INS results have been interpreted [0 in 
terms of the additional bound state contribution predicted in Ref . Thus, our results may 
not be of direct relevance to this particular experimental system. However, our work does 
predict that a second peak in the INS cross-section should be seen in strongly coupled ladder 
systems and it is quite possible that the feature persists to all orders in the perturbation 
expansion we have employed. It would be interesting to confirm this effect by looking at 
other systems that are more convincingly modelled by a simple ladder Hamiltonian and it 
is hoped that future experiments do indeed see the effects coming from the bound state. 



C. Weakly-coupled two-leg ladders 



In this section, we analyze the ladder system fl2.4| ) in the complementary weak-coupling 
limit: J± <C J\\. An elegant mapping developed by Shelton et.al. |25|] allows one to express 



the low-energy, long-distance properties of the model in terms of a continuum theory of 
weakly-interacting massive Majorana (real) fermions. We will analyze the low-energy scat- 
tering properties of the spin one excitations of the weakly-coupled ladder by working in this 
Majorana fermion representation. 

We begin with a brief review of the Majorana fermion representation. We will not at- 



tempt here to describe in any detail the procedure used |25j to arrive at this field-theoretic 
representation. Instead, we will be content with a rather telegraphic summary of the prin- 
cipal steps involved. To begin with, one writes down the usual, free, massless bosonic 



theory [13] for the low-energy properties of each of the two S = 1/2 Heisenberg antiferro- 



magnetic chains that make up the ladder. The interchain exchange J± is then turned on, 
introducing a local, isotropic (in spin space) coupling between the spin-density operators 
of each chain in the bosonic representation. This has two pieces to it: one coupling the 
staggered parts of the spin densities with each other and the other doing the same for the 
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uniform component. Now, one works with symmetric and antisymmetric combinations of 
the two boson fields (one for each chain) and transcribes everything to a fermionic repre- 
sentation, introducing one Dirac fermion for the symmetric combination and another for 
the antisymmetric combination in the usual manner (for a readable account of the relevant 
machinery of Abelian bosonization, see for instance the review |RJ by Shankar). The last 
step is to write each Dirac fermion as two Majorana fermions. If one leaves out the uni- 
form part of the coupling to begin with, the theory in terms of the Majorana fermions is, 
remarkably enough, a free-field theory. The staggered part of the coupling just provides a 
mass A to each of the two Majorana fermions obtained from the symmetric combination of 
the bosons, while the two Majorana fermions obtained from the antisymmetric combination 
acquire masses A and— 3A respectively (the actual energy gap is given by the absolute value 
of the mass). The three Majorana fermions with mass A form the spin one triplet we expect 
on general grounds, and the fourth Majorana fermion represents a high-energy singlet mode 
that will not be very important for our purposes. The mass parameter A of the theory is 
proportional to J± with the proportionality constant being non-universal. Finally, turning 
on the coupling between the uniform part of the spin densities gives us a four-fermion inter- 
action term between these massive Majorana fermions which will play a crucial role in our 
analysis of the iS-matrix. 

The procedure outlined above gives us the following Hamiltonian for the effective field 
theory written in terms of Majorana fermions: 

H= £ H A (n+n^ A (p)+Hr; (2.22) 

a=x,y,z 

here the £ a and p are Majorana fermion fields with anticommutation relations given as 

{pR(x),Pn(y)} = S(x-y) , 
{p L (x),p L (y)} = 6(x-y) , 
m^),e R (y)} = S ab 5(x-y), 

{Q{x),&(y)} = 5 ab 6{x-y) , (2.23) 
with all other anticommutators being equal to zero, 7^ m (0) is defined in general as 

Wm(0) = -y- J dx{(j) L d x (j) L - 4> R d x (f> R - m(p R (f) L ) , (2.24) 
with vf ~ J\\a and the interaction term Tii may be written as 

n^gj dx{c R ae R a + aae R e L + mtiti - mi + e R a + w*} , (2.25) 

with g ~ Jj_a. Note that each Majorana fermion is a two component object, the two 
components being labelled with the subscripts R and L to denote the 'right-moving' and 'left- 
moving' parts. To make contact with the original spin-ladder, we also need a prescription 
for expressing the spin-operators of the ladder in terms of the Majorana fermions. In sharp 
contrast to the a-model, only the uniform part of the spin-density operator has a local 
representation in terms of the fermions; the components of the spin-density near q = Q can 
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be expressed only in terms of highly non- local functions of the fermi- fields ||25|| . We have the 



following expressions |25[ for the uniform parts, J\ and J2, of the spin density on each chain: 

i ,1 



i ,1 



(*) = ^ abc C{x)C{x) - C{^)Pu{x)) , (2.26) 

where the index v takes on values R or L and repeated indices are summed over. Note that 
the field p corresponding to the non-universal high-energy singlet mode drops out of the 
expression for the uniform part of the total spin-density of the ladder which can then be 
expressed entirely in terms of the spin one triplet fields. 

We shall find it convenient, when it comes to actually doing any calculations, to rewrite 
all of the foregoing in terms of fermionic creation and annihilation operators. These are 
defined as follows: Let and p v {p) denote the Fourier transforms of and p v (x) 

respectively. We write 

Pu(p) = 9u{p)s{p) + g~v{-p)s\-p) , (2.27) 

where t a (p) and s(p) are the fermionic annihilation operators corresponding to the triplet and 
singlet modes respectively and f v (p) and g u {p) are complex- valued functions of p which we 
specify below. These creation and annihilation operators obey the usual anticommutation 
relations: 

{t a (p),tl(q)} = 2n5 ab 5(p-q) , 

{s(p),s\q)} = 2n5(p-q) , (2.28) 

with all other anticommutators equal to zero. In terms of these operators, the non-interacting 
part of the Hamiltonian reads 

n = r %(p)ti{p)t a {p) + r % s {p)s\p) S {p) , (2.29) 

J-00 Z7T J-00 /7T 

where e(p) = (p 2 v 2 F + A 2 ) 1//2 , e s (p) = (p 2 v 2 F + 9A 2 ) 1 / 2 , and the repeated index a is summed 
over. The functions f u and g v are actually chosen to ensure that the non-interacting Hamil- 
tonian has this simple diagonal form in terms of the creation and annihilation operators; 
this choice guarantees that the operators and t\, as defined in (|2.27|) , create the true 



quasiparticles of the non-interacting system. The expressions for f u and g v are best written 
as follows: 

!r{p) = «a(p) P > , 
fn(p) = iv A (p) p < , 
hip) = fn(-p) Vp , 

9 R^) = «(-3A)(p) P>0, 

9r(p) = ««(-3A)(p) P<0, 

9L(p)=g R (-p) Vp; (2.30) 



18 



here the functions u m (p) and v m (p) are defined in general as 

u m (p) = cos(0 m (p)/2) , 

v m (p) = sin(0 m (p)/2) , (2.31) 

with the angle 9 m (p) being specified by cos(9 m (p)) = VF\p\/(m 2 + vpp 2 ) 1 ^ 2 , sm(6 m (p)) = 
msgn(p)/(m 2 + Vpp 2 ) 1 ^ 2 . Now, we can rewrite the interaction term in normal ordered form 
with respect to these singlet and triplet creation and annihilation operators. The quadratic 
terms so generated give the first order correction to the masses of the singlet and triplet 
modes (this correction has already been calculated in Ref [25| by other means). The quartic 



term left over, has, in addition to the usual, normal-ordered, particle-number conserving 
piece, other pieces that involve pair creation and destruction. The full expressions are 
somewhat messy and we refrain from displaying them here. However, and this is key, 
we will need only a very simple part (corresponding to the low momentum limit of the 
particle-number conserving piece) of this quartic term for the calculation of the iS-matrix 
in the low- momentum limit; our method of writing everything in terms of the creation and 
annihilation operators has the advantage of identifying and isolating this piece at the very 
outset. Finally, as an aside, we note that the total spin operator of the system may be 
written in terms of the triplet operators as 

Sit = ^ bC ^tt(p)t c (p) ; (2.32) 

this confirms that the triplet creation operator t\ does indeed create a single spin one quasi- 
particle (with polarization a) of the non-interacting system. 

With all of this in mind, let us turn to the analysis of the scattering properties of this 
model. As we are hoping to calculate the iS-matrix perturbatively in the coupling g, it 
is convenient to write 5 = 1 + iT . The 'transition-matrix' T can then be calculated 
perturbatively using the standard field-theoretic prescription that relates it to the corre- 
sponding amputated, connected Green's functions of the theory. Let us make this precise 
for the case we are interested in: namely, a scattering process in which the initial state 
consists of two particles, one with momentum k\ and spin polarization and the other 
with momentum k 2 and spin polarization /i 2 , and the final state has two particles labelled 
by (k[, [i'i) and {k' 2 ,n 2 ). Note that we are now not talking about the bare particles of the 
non-interacting theory, but the actual physical quasiparticle states of the system, correct to 
the relevant order in the perturbative expansion in g. The corresponding matrix element, 
S£fff(ki, k 2 ; k[, k' 2 ) = (k[ k' 2 n' 2 \S\ki fj,i, k 2 /x 2 ), may then be written as 

Sggfa, h; k[, k' 2 ) = (2 7 r)% lK 5^ 2 5(k 1 - k[)5(k 2 - k' 2 ) + 

(2n) 2 6(E f - Ei)S(k f - k^iM^ih, k 2 ; k[, k' 2 ) , (2.33) 

where Ef = e(k[)+e(k' 2 ) and = e(ki)+e(k 2 ) are the final and initial energies respectively, 
kf and hi are the total momenta in the final and initial states respectively, and A4 is the 
'reduced' matrix element (with energy and momentum conserving S functions removed) for 
the process under consideration. 
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We now specialize to the case k\ = k, k 2 = —k (k > 0); this special case allows us to make 
our basic point (regarding the infrared divergences present in a perturbative calculation 
of the scattering properties) while keeping the calculations simple. In this case, we may 
decompose the scattering matrix as follows: 



s£2(*> K, K) = 5(k[ - k)5(k' 2 + k) \s x {k)5^8 t 



S2(A0<W<W + S 3 (k)S^5^ . (2.34) 



Now, energy and momentum conservation in one dimension provide enough constraints on 
the two-body problem to ensure that the allowed final states have the same set of momentum 
labels as the initial state. This allows us to convert the overall energy and momentum 
conserving S functions in the second term of (|2.33|) to S functions that identify k[ with k 



and k! 2 with —k. In the process, we of course introduce additional kinematic factors coming 
from the Jacobian (we are basically using 5(f(x)) = 5(x)/\f'(x)\ ). Using this, we can write 

*<*) = ■ 



E(k) 
>F 

where Mi, M 2 , and M 3 are defined in terms of the following decomposition for M. 



Ss(k) = ) iM 3 (k) , (2.35) 



K\7'S k ' ~ k > k > = KiWWA^ + ^(A0<Wi<W + ^(*) . (2.36) 



The relations ( |2.35| ) are useful because there is a simple diagrammatic prescription for the 



perturbative evaluation of M.. According to this standard field theoretic prescription p6 
i-M-^^ki, k 2 ; k[, k' 2 ) is proportional to the sum of all 'amputated' (factors corresponding to 
external legs omitted), fully connected, one particle irreducible diagrams contributing to the 
time ordered four-point function with two incoming external lines and two outgoing external 
lines. The incoming lines must carry momenta ki and k 2 , frequencies u)i and uj 2 set to their 
respective 'on-shell' values of e(ki) and e(k 2 ), and spin labels /ii and fi 2 respectively. The 
outgoing lines must carry momenta k[ and k' 2 , frequencies again set to their on-shell values 
of e(k[) and e(k' 2 ), and spin labels fi[ and ji' 2 respectively. Denoting the sum of all such 
diagrams schematically by 1^, we can write 

iM Zk^ k ^ = (v^) 4 r4(fci /ii, k 2 // 2 ; k[ fi[, k' 2 //' 2 ) , (2.37) 

where the field-strength renormalization factor Z comes into play because the singlet and 
triplet creation operators and t\ create the bare particles, while we are asking questions 
about the scattering properties of the physical quasiparticle excitations. We will not be very 
careful here about the precise definition of Z; it will soon become apparent that this does 
not play any role in the calculation we do. 

Before we set about calculating T 4 , we need to specify our conventions regarding the 
diagrammatic representation of perturbation theory. As shown in Fig |5], we denote the 
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propogator of the triplet particle by a solid line with an arrow carrying momentum k, 
frequency u and spin label //; this has a factor of if (u — e(k) + irj) associated with it. It 
turns out that we do not need to consider any diagrams that have lines corresponding to 
singlet particles and we will not bother to introduce a diagrammatic representation for their 
propogator. We also display our diagram convention for the four point vertices of the theory 
in the same figure; again, only the particle number conserving vertices in which all four 
lines involved correspond to triplet particles have been assigned a diagram as the others 
will not play a role in what follows. One type of vertex, labelled (a) in Fig |^, depicts a 
process in which two particles of momentum p 3 and p±, both with spin label fi = x scatter 
into a final state populated by two particles with momenta p\ and P2, and spin label fi = y. 
The full momentum dependent factor associated with this diagram is also shown below it. 
We will need only a very simple low-momentum limit of this expression in most of what 
follows. The other kind of vertex, labelled (b) in Fig |5], shows incoming particles with labels 
(p 3 y) and (p 4 x) scattering into a final state populated by particles with labels (pi x) and 
[j>2 y) respectively. Again, the full momentum dependent factor is displayed alongside for 
completeness. We will mostly need only the value of this factor when all four momenta 
equal zero; this is given simply by —ig. Of course, all other vertices of the same type, but 
having different spin labels that can be obtained from these using the 0(3) symmetry of the 
problem, have the same factors associated with them. 

We are now in a postition to do some calculations. We begin by noting that, apart from 
the overall factor of (v^) 4 which we are ignoring for now, iM\{k), iM 2 (k) and iM 3 (k) are 
equal to r 4 (/c x, —k x; k y, —k y) , r 4 (/c x, —k y; k x, —k y) , and r 4 (&; x, —k y; k y, —k x) re- 
spectively. It is quite simple to calculate these three quantities to leading order in g. The 
diagrams contributing to zMiis shown in Fig |6], while those contributing to iM 2 and iM 3 are 
shown in Fig [7[ Evaluating these 'tree-level' amplitudes, we obtain 

k 2 v 2 

T 4 (k x, -k x; k y, -k y) = W^nr > 

A 2 

T 4 (A; x, -k y; k x, -k y) = W^JT , 

T 4 (fc x, —k y\ k y, —k x) = —ig . (2.38) 

As long as we are interested in only the first order result for S, we can set Z — 1 and directly 
use these expressions to get the following results for the leading low k behaviour of Si, 5*2, 
and S3 correct to first order in g: 



Si(k) 



ig ( kvF^ 



S 2 (k) = 1 + 



2v F v a ; 
ig ( A 



2vf \kvF 



2vf \kvF 

We immediately see that the perturbative expansion cannot be trusted in the low-momentum 
limit because of the infrared divergences present in the expressions for 5*2 and 5*3. The 
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structure of this first order result is seen to be qualitatively similar to the 0(1/N) result for 
the two-particle iS-matrix of the O(N) cx-model [fLOH . In the latter case, we know that the 
exact value of the iS-matrix is perfectly well-behaved in the k — > limit and is in fact given 
by the super-universal expression (|1.2j ). To obtain the correct result in this limit for our 
problem, we need to identify the leading infrared divergences at each order in g and perform a 
resummation. Now, we do not expect any infrared divergences in the perturbation expansion 
of Z and as a result the pref actor of (VZ) 4 in the expression for A4 does not contribute at 
all to the terms that need to be resummed; we will forget about this factor from now on. 

Let us now try and identify the leading infrared divergent diagrams at each order in 
perturbation theory. First of all, it is clear, purely from frequency and momentum conserva- 
tion at each vertex, that no diagrams involving pair creation or annihilation can provide the 
leading divergence at any order. Moreover, only internal loops in which both propogators 
involved point in the same direction give a nonzero result on doing the integral over the 
frequency running through the loop. A little thought should convince the reader that these 
two constraints allow us to conclude that the ladder series shown in Fig || give the leading 
infrared divergent terms in 5*2 and S3 to all orders in g. Turning our attention to Si, we 
see immediately that Fermi statistics guarantees that each vertex in the analogous ladder 
series for Si has enough factors of momentum associated with it to rule out any infrared 
divergence appearing in Si. Our task is thus reduced to evaluating the two series shown in 
Fig H To do this, we note that as far as the coefficient of the divergent piece is concerned, 
we can ignore the momentum dependence of each vertex and simply replace it with a factor 
of — ig. Each crossing of the fermion propogators gives a factor of —1 and each loop integral 
gives A/2kv F . Putting all this together and summing the resultant geometric series, we 
obtain the following non-perturbative results for the low momentum behaviour of S 2 and S 3 : 

s 2 (k) - 2lkv2p 



gA + 2ikv F 
gA + 2ikv' F 



Ss(k) = • (2.40) 



An interesting feature of these results is the pole in the upper-half k plane at k = iAg/2v F 
which seems to suggest the presence of a bound state. However, this region of k space is 
definitely beyond the domain of validity of (|2.40|) and it is not clear what significance, if 
any, to ascribe to this curious fact. 

Turning to firmer ground, we see that the foregoing implies that the low-momentum 
limit of the two particle 5-matrix is perfectly well-defined and is in fact given by 

S^(k h k 2 ; k[, k' 2 ) = S^ 2 S^2nS(h - k[)2n5{k 2 - k' 2 ) . (2.41) 

Note that apart from an overall factor of minus one, this is exactly the super- universal form 
( |1 . 2|) . The relative sign is simply a consequence of fermi statistics and our choice of phase for 
the final state of the scattering process. In any case, we will see that when we use the super- 
universal form of the S- matrix for discussing spin transport, the overall phase is immaterial. 
On the other hand, the overall factor of —1 in the superuniversal form ( |1.2| ) will be crucial 
when we work out the correlators of the staggered component of the magnetization density. 
This may seem worrisome at first sight. However, as we do not have any local representation 
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of the staggered component of the spin density in terms of the Majorana fermion operators, 



there is no contradiction at all. In fact, the semiclassical techniques used in Section [[II A 
may also be applied to the problem of calculating the finite temperature Green function of 
the fermions; this would correspond to calculating the finite temperature correlators of some 
highly non-local string operators of the original spin system. However, as it is difficult to see 
how these may be accessible at all to any experimental probes, we do not pursue this line 
of thought any further. 

Thus, we see that the low-momentum behaviour of the 5-matrix in this fermionic rep- 
resentation of the weakly-coupled ladder is consistent with the super-universal form ( |1 . 2| ) , 
although this behaviour is definitely not accessible to perturbation theory. This leads us 
to believe that similar infrared divergences would invalidate any perturbative calculation of 
dynamical properties at finite temperature (when there will be a dilute gas of thermally ex- 
cited particles present) that uses this representation. In particular, this appears to indicate 
that the results of Ref [EH] for the NMR relaxation rate \jT\ are incorrect at low T > 0. 



III. DYNAMICS AND TRANSPORT FOR 0<T<A 

The results of this section are expected to apply to all gapped one-dimensional antiferro- 
magnets with massive spin-one quasiparticles. We will develop, what we believe is an exact 
semiclassical theory of dynamics and transport for T < A. We will consider fluctuations 
near q = Q in Section III A, and near q = in Section III B. 



A. Thermal broadening of the single-particle peak in S(q,tv) 

In this section, we present calculations leading up to our results for the thermal broad- 
ening of the single particle peak in the dynamics structure factor. 

The inelastic neutron scattering cross-section provides a direct measure [jl| of the dy- 
namical structure factor S(q, uj) which is defined as 

S(q,u) = ^J dte^(s a (q,t)s a (-q,0)}; (3.1) 

where s a (q,t) is the Heisenberg representation operator corresponding to the component 
of the spin density at wavevector q, the expectation values are with respect to the usual 
equilibrium density matrix and summation over the repeated index a is implied (note that 
we are assuming rotational invariance in spin space and working at H = 0). We are interested 
in the structure factor for q close to ir/a. In this case we have 

S(q,u) oc — / dte iujt (n a (k,t)n a (-k,0)); (3.2) 
2tc J 



where k = q — ir /a. To get a feel for what ( |3.2|) looks like at T = 0, let us consider a 



particular lattice regularization of the cr-model, defined by the quantum rotor Hamiltonian 

z i 9 i 
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where Lj is the angular momentum operator of the rotor at site i, hi is the unit vector that 
denotes the position of the rotor at site i and we have temporarily set c = a = 1. It is 
not hard to analyze the properties of this model in a large g, strong coupling expansion; 



moreover this is expected to be qualitatively correct for all g in one dimension ||47|| . To 
lowest order, we can easily see that the ground state would just be a product state with 
each site being in an eigenstate of L with zero eigenvalue. The lowest excited states would 
be a degenerate manifold corresponding to promoting any one site to the L = 1 state and 
thereby creating a 'particle' at that site. To first order in 1/g, a hopping term would be 
generated in the effective Hamiltonian for the single-particle sector, resulting in a band of 
one particle excitations. To this order, n is just a sum of creation and annihilation operators 
for the stable particle-like excitation of the system. At higher orders in 1/g, n acting on the 
vacuum will also produce multiparticle states, but there will always be some single particle 
component. Reverting back to our continuum theory, we see that ( |3.2| ) evaluated at T = 
would have a contribution ~ 5{uj — s(k)) associated with the stable particle. The next 



contribution is actually a continuum above the 3-particle threshold fUj . Following |23], we 
shall now focus exclusively on how this one-particle peak broadens as T becomes non-zero. 
Let us define 

C(x,t) = -(n a (x,t)h a (0,0)) , (3.3) 

where the repeated index a is summed over. Let K(x,t) denote C(x,t) evaluated at T = 
keeping only the single particle contributions. We have 

K(x,t) = I ^D(p)e ipx - i£{p)t . (3.4) 

J Z7T 

Here D(p) is a 'form factor'. For our Lorentz invariant continuum model, 

D ^ = W) (3 ' 8) 

where A is a non-universal quasiparticle residue. This gives K(x,t) = AK (A(x 2 — 
cH 2 ) 1 / 2 1 c) I (2ir) , with K the modified Bessel function. 

Now let us evaluate C (x, t) for non-zero temperatures using the semiclassical method of 
2"3fl . First, it is convenient to switch to operators n + i(x), n_i(x) and tiq(x), defined as 



n+i = nli = n x - in y , 



and 



n = n z . 

n+i is a sum of a creation operator for particles with z-component of spin m equal to +1 
and an annihilation operator for particles with m equal to —1. no is a sum of creation and 
annihilation operators for particles with m equal to 0. In the absence of an external field, 
we may write 

C{x,t) = (n {x,t)no{0,0)) . (3.6) 
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We represent ( |3.6|) as a 'double time' path integral, with the e lHt factor coming from the 
Heisenberg operator generating paths that move forward in time, and the e +lHt producing 
paths that move backward in time. We begin with an initial state which is populated by 
thermally excited particles, the density of particles being ~ e~ A/,T and their mean spacing 
being much larger than the thermal de-Broglie wavelength ~ c/(AT) -1 / 2 . As argued in 
22].E3|, this means that the particles can be treated semiclassically. In this semiclassical 



limit the dominant contribution to the Feynman sum comes about when the paths going 
backward in time are exactly the time-reversed counterparts of those going forward and all 
particles follow their classical trajectories between collisions p2"| , |23|| . Whenever two particles 



collide, energy and momentum conservation is sufficient to determine the final momenta. 
However, one cannot entirely ignore quantum effects of the collisions. The spins of the 
particles after the collision as well as the phase picked up by the wavefunction of the system 
as a result of the collision is determined by the quantum mechanical scattering matrix (S). 
For T < A, the particles all move very slowly and we need only the super-universal low- 
momentum limit (O). 



All this leads to the following description of C(x, t) in this asymptotic limit [22]: At time 
t = we begin with an initial state populated equally (for H = 0) with three species (cor- 
responding to the three values of spin projection m) of particles each uniformly distributed 
in space with density p/3, where the total density p is given as 



P 



3 f ^Pg-(A+cV/2A)/T = 3 / TA -A/T 

J 2tt V 2ttc 2 



The velocities are distributed according to the classical Maxwell-Boltzmann distribution 
function 



Each particle in the initial state is assigned one of the three values of m with equal probability 
(assuming H = 0). The operator no(0) acting on this initial state creates at time t — 
one extra particle at x = with spin value equal to (the annihilation part of no gives a 
contribution which is exponentially suppressed and is ignored here). These particles follow 
their classical trajectories forward in time. At every collision, we pick up a factor of — 1 
from the iS-matrix. At time t, a particle with spin projection of zero is annihilated at x by 
n (x). The resulting state is then propogated backward in time to t = and its overlap 
with the initial state calculated. C(x, t) is then given by the average of this overlap over the 
ensemble specified earlier. 

A typical example of a space-time configuration of trajectories that leads to a non-zero 
value for this overlap is shown in Fig |9]. All trajectories in the figure except the dotted line 
denote space-time paths that are traversed both forward and backward in time. The dotted 
line is traversed only forward in time as the particle travelling on it is destroyed at time t by 
no(x). A little thought convinces one that this overlap is non-zero only when all particles 
colliding with a particle travelling on the dotted trajectory have the same spin m (equal to 
zero) as it does. Moreover, when this condition is satisfied, the value of the overlap is just 
(— l) Ul K(x, t) where rii is the number of collisions that the dotted trajectory suffers. The 
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factor of (— l) n ' comes from the scattering matrix at each collision between a particle on the 
dotted trajectory and other particles. All other collisions occur in pairs (the second member 
of the pair coming from the evolution backward in time) and thus do not contribute any 
phase factor. The factor of K(x,t) is just the relativistic amplitude for the propogation of 
a single particle from x = at t = to position x at time t. 
All this implies that we can write 

C(x,t) = R(x,t)K(x,t) (3.7) 

which defines the 'relaxation function' R. For the case where the particle has only one 
allowed value of its spin label, m, it is possible to compute R(x,t) analytically ||23|| , and the 
resulting expression ( |3. 7] ) then agrees precisely with a computation using very sophisticated 
quantum inverse scattering method |27j : this agreement gives us confidence that the physical 



approach developed here is asymptotically exact at low temperatures. 

Let us now turn to the calculation of R for the case of interest here. We begin by writing 
a formal expression for R based on the foregoing semiclassical description. Let {x k (0)} be 
the positions of the thermally excited particles at time t — 0. Let {v k } be their initial 
velocities. Here k is an index running from 1 to N, the total number of particles present 
in the initial state in a system of size L. We label the initial positions with the convention 
that x k (0) < xi(0) for k < I. Let X k (t) = x k (0) + v k t denote the k th space-time trajectory 
(note that this is quite distinct from the position of the k th particle at time t). Let m k (t) 
denote the spin projection value of the particle travelling along the k th trajectory at time 
t. The spin projections are randomly assigned to each particle at time t = as described 
earlier and m k (t) at later times depends on which particle is travelling on the k th trajectory 
at any given time. We have the following expression for R: 

R(x,t) = ([[r k ) , (3.8) 

k 

with 

Fk — 1 — (1 + 5m k (T k ),i)®k ; 

where 

6 fc = 6{x - X k (t))d(x k (0)) + 6(X k (t) - x )6(-x k (0)) , 

and 

r k = x k (0)t/(x - v k t) . 



The angular brackets in ( |3.8| ) denote averaging over the ensemble of initial conditions spec- 
ified earlier. 

Unfortunately, it does not seem possible to do the ensemble average analytically. Using 



methods of Refs [f49| , [50ll i it is possible to develop a 'cumulant' expansion for the logarithm 



of R |[51|| . This expansion, however, is essentially a short time expansion which is not 
uniformly convergent, and thus not very useful for our purposes as we eventually need to 
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Fourier transform C(x,t). It is also possible to develop a 'mean-field' approximation to 
this classical model that ignores the complicated correlations between the m&(t) at different 
times (see Appendix D). This proves to be reasonably accurate at least for the R(0,t), 
though the high-accuracy numerics we describe next show clear deviations from the mean- 
field results. So, although we have an asymptotically exact formulation for the non-zero 
temperature C(x,t) at distances much larger than the thermal de-Broglie wavelength and 
times much larger than T -1 , we need to numerically determine the relaxation function R to 
actually calculate anything accurately. This is what we turn to next. 

An important property of R(x,t), which follows directly from fl3.8| ) is that it can be 
written in a scaling form as R(x, t) = R(x, t) with x = x/L x and i = t/L t where 

r -I 



P 



and 



1/ a y/ 2 

\2c 2 TJ 



P 



Thus it is most convenient for the numerics to measure length in units of L x and time in 
units of L t and directly calculate R. We start with a system size of L — 400 (in units 
of L x ) and impose periodic boundary conditions. The density in these units is unity and 
so the initial state is populated by 400 particles with their initial positions drawn from a 
uniform ensemble. This system size is large enough that finite-size effects are negligible for 
our purposes. Each particle is assigned a velocity from the classical thermal ensemble. In 
these new units this implies that we choose velocities from the distribution 

V(v) = ±=e-* 2 . 

'71 



An important advantage of our method is that we do the average over the spin values 
analytically. To do this, we note that it is possible to reformulate the calculation of R by 
writing 

R = (T(C)> , 

where C denotes a given space time configuration of trajectories, the angular brackets denote 
averages only over the initial positions and velocities that define this configuration, and T(C) 
is defined as 



T(C) = (-!)•* Q' 



Here, nh is the total number of collisions involving a particle travelling on the dotted tra- 
jectory of Fig y and n p is the number of different thermally excited particles that have 
had collisions with a particle travelling on the dotted trajectory. Now, T(C) = for all 
configurations C in which the presence of the extra particle (that starts out on the dotted 
trajectory) affects the evolution of the various spin values rrikit). So we might as well forget 
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about the particle travelling on the dotted trajectory and consider an auxiliary spacetime 
diagram that now involves only the thermally excited particles. We now agree to ignore the 
spin label on the dotted line of Fig |9|; the dotted line now does not denote the trajectory of 
any particle. In terms of this picture we can define as the number of times any solid line 
crosses the dotted line, and n p as the number of different thermally excited particles that 
cross the dotted line. 

With this new formulation in hand, calculating T(C) reduces to some simple book- 
keeping that keeps track of these two integers for a given configuration C. We implement 
the ensemble average by averaging over 4 x 10 6 configurations drawn from the appropriate 
distribution. The combined absolute error (statistical and finite-size) in R(x, t) for values of 
x, t of interest to us is estimated to be less than about 5 x 10~ 4 . 

With R available, it is a simple matter to numerically Fourier transform the resulting 
C(x,t) and obtain the dynamic structure factor S(q,u). Details of the numerical procedure 
used are relegated to Appendix B. Here we only comment on some conceptual issues involved 
and discuss our results. 

There is an important subtlety associated with doing the Fourier transform that needs 



to be first addressed. As discussed in Ref [p3| , the semiclassical result for C(x,t) is valid as 



long as both x and t are not very small; the results break down when x ~ Xt and t ~ 1/T 
(At being the thermal de-Broglie wavelength). However, the Fourier transform of C (at 
wavevector k = q — n/a) is an asymptotically valid approximation to S(q,u) only for u close 
to e(k). The reason for this can be understood by noting that the long-time asymptotics of 
our form for C(x,t) have an oscillatory character with oscillations on the scale of A -1 . Put 
another way, it is the spectral weight in the one-particle peak that plays a dominant role in 
determining the long-time, large- distance asymptotics of C(x,t) and so we can learn only 
about this feature in the spectral weight by Fourier transforming our form for C . 
With this caveat in mind, we have 

S(q,u) = — I I dtdxe iujt - ikx K(x,t)R(x,t) . (3.9) 

where k = q — it /a. We have not attempted to exhaustively map out S(q,ui), although 
it would be quite straightforward to get more extensive numerical results should they be 
of interest in some experimental context. Below we confine ourselves to discussing our 



results for S(q,tu) for a couple of sample values of q. Fig |10| shows scans in frequency 
across the quasiparticle peak in S(q,u) forq = Q at four different values of temperature. 
It is interesting to note that when rescaled by L t and plotted against a rescaled frequency 
variable 5lj = (uj — A)L t , the three curves for A/T = 3, 4, and 5 seem to collapse on top 
of one another within our numerical errors (which are conservatively estimated to be a few 



percent at the most). In Fig [11], we show a scan in wavevector across the same peak for 
lo — A, again at the same four values of temperature. The curves at the lower temperatures 
again show scaling collapse; when rescaled by L t and plotted against the rescaled variable 



k = kcyL t /A, they seem to all fall on top of one another. Moreover, the scaling curve 

in Fig JTj, when plotted as a function of the independent variable —k 2 /2 coincides within 
our numerical error with the scaling function of Fig |10| for 8u < 0; this is displayed in 



Fig O. While we do not have any reason to expect that this scaling is generally true, all 
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three observations may be put together in terms of a scaling form that is valid locally in the 
neighbourhood of the quasiparticle peak for q = Q; more formally we write 



AcL t (lo — e(kY 
~^A I Lt 



. (3.10) 



We also investigated S(q,co) in the vicinity of the quasiparticle peak corresponding to 
q = Q + A/c ; for this to be meaningful, we of course need A/c to be much less than the 
microscopic scale ~ a -1 beyond which our continuum theory does not work. We again tried 
to check if the analogous scaling form, 

„, . AcL t ^ ( uj — e(k)\ /„„„v 



is approximately valid. Fig [13] shows scans in frequency across the peak with k held fixed 
at A/c, for A/T = 2, 3, 4, and 5. We see that the curves do not really appear to scale. 
In Fig [14], we show scans in wavevector, with lo held fixed equal to v^2A for the same 
values of the ratio A/T. We plot the data (rescaled by L t ) against the rescaled variable 
5k = cL t (k — A/c) (note the difference in the choice of rescaling of the independent variable 
from the earlier case). Again, in sharp contrast to the q = n/a peak, we see that the curves 
do not show any signs of scaling; our local scaling form is not a very good way of organizing 
the data in this case. 

These scaling properties are best understood as follows: Imagine developing R(x, t) in an 
expansion about x = for constant t and then calculating the Fourier integral in (|3.9| ). The 
zeroth order term clearly gives us a result for S(q,u>) which is compatible with the scaling 
form we have postulated for asymptotically low temperatures. However, before we can trust 
this result, we need to check that the corrections to the leading behaviour go to zero in the 
limit T — > 0. This is where the difference between the two peaks we looked at becomes 
apparent. It is easy to see that this is true only for values of q such that c 2 \k\L t /AL x — > 
as T — * and this explains why the scaling form ( |3.11|) does not work. Now consider 



the peak at q = Q: The zeroth order scaling result has most of its weight in the region 
\k\ < \/~A/c\f77 t . For \k\ ~ vA/cy/Lt, the corrections to this leading result do indeed go to 
zero and this establishes the scaling form fl3.i0|) . An interesting feature of this result is that 
the scaling function $ is completely determined by the x = part R(0, t) of the relaxation 
function. More precisely, we have 

$(z) = - f° dse izs R(0,s) . (3.12) 

A useful check on all of our calculations is thus to compare the scaling function obtained 
in Fig [T(] and Fig [12] with ( |3.12| ) evaluated numerically (it is possible to do this to a high 



accuracy; details may be found in Appendix B). The results of such a comparison are shown 
in Fig 13 and the agreement is seen to be quite good. While the numerical results for R(0, t) 



show a clear deviation from the simple exponential decay predicted by the 'mean-field' theory 
referred to earlier, we do find that the corresponding simple Lorentzian form for the Fourier 
transform: $(z) = 7ta/2(a 2 + z 2 ) (with a ~ 0.71) provides an excellent approximation to 
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the line-shape (the 'mean-field' theory, however, gives a value of 4:/3^/^T m 0.7523 for a — see 
Appendix D). 

We thus have results for the thermally broadened quasiparticle peak in S(q,u>); the 
accuracy of these in the asymptotic regime (T -C A) is limited only by the computer time 
spent in numerically evaluating the relaxation function and doing the Fourier transform. 
These results, especially the scaling properties in the vicinity of the peak corresponding to 
q = Q, should be of relevance to neutron scattering experiments on gapped one-dimensional 
Heisenberg antiferromagnets performed at temperatures T <C A and it is hoped that this 
study provides a useful paradigm for organizing the experimental results. 



B. Low temperature spin diffusion probed by 1/T\ 

In this section, we shall present a detailed comparison of our results J22| for the field 
(H) and temperature (T) dependence of the NMR relaxation rate 1/Ti (in the regime 



T, H <C A) with the experimental data of Ref [3(| on the NMR relaxation rate in the 



compound AgVP2S6 which is thought to be a S = 1 one-dimensional antiferromagnet with 



a large gap A m 300 K and single-ion anisotropy energy of about 4.5 K ||30|| . We will 
ignore this anisotropy for the most part in our theoretical analysis (although we are forced 
to phenomenologically introduce spin-dissipation into our theorerical results in order to fit 
the data of Ref |30| at low temperatures, we do not have any theory that gives the detailed 



temperature dependence of this spin dissipation rate starting from the anisotropic coupling 
term in the Hamiltonian) . 

For completeness, let us begin with a detailed review of the calculations leading up to 
our expression for 1/Ti. The NMR relaxation rate is given in general by an expression of 
the form 

1 r dq 

7f7=Y, £ / ^Z A ap{(l)A ai (-q)Sp 1 (q 1 uj N ) ; (3.13) 

J l a=x,y f3,j=x,y,z J Z71 ~ 

where S^(q, u) is the Fourier transform of the spin-spin correlation function (the subscripts 
refer to the 0(3) indices of the spin operators), uj^ = JnH is the nuclear Larmor frequency 
(7at is the nuclear gyromagnetic ratio), the field H points in the z direction and A a p are 
the hyperfine coupling constants. The q integral in ( 3.1 3|) is dominated by values of q near 



|U) and we can thus work out the field and temperature dependence of 1/Ti knowing the 



T > correlators of the conserved magnetization density of the 0(3) NLcrM field theory. 
This is what we turn to next. 

We define the correlation functions 

C u>zz {x,t) = (L z {x,t)L z {0,0)) - (L) 2 

C u ^ + (x,t) = (L_(x,t)L + (0,0)> ; (3.14) 

here the angular brackets denote averaging over the usual equilibrium density matrix, L z (x, t) 
is the Heisenberg representation operator corresponding to the z component of the magne- 
tization density, and L± are operators corresponding to the circularly polarized components 
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of the magnetization density defined as L± = L x ±iL y . As argued in Refs [£^,[23|], these cor- 
relation functions in the asymptotic regime may be evaluated by writing down a double-time 
path integral representation for them and evaluating it semiclassically. 



This leads to the following prescription [2~2~ for C UiZZ (x,t) : At time t = we begin 



with an initial state populated with three species (corresponding to the three values of 
spin projection m) of particles each uniformly distributed in space with densities given 
respectively by 



dP_ -(A-mH+cV/2A)/T _ M^.„-(A-mH)/T 

2tt V 27T c 2 



and with velocities distributed according to the classical Maxwell-Boltzmann distribution 
function 



Each particle in the initial state is assigned one of the three values of m with probability 
fm — e mH l T /(I + 2cosh(i?/T)) . The operator L z (0) merely keeps track of the local value 
of the z component of the spin. Acting on the initial state, it measures the z component of 
the magnetization density in the initial state at position x = 0. These particles then follow 
their classical trajectories forward in time. At every collision, the particles retain their spin 
labels. In addition, the state picks up a factor of —1 from the S- matrix at each collision. At 
time t, the operator L z (x) measures the value of the z component of magnetization density 
at position x. The state is then propogated backward in time to t = and its overlap with 
the initial state calculated. C UjZZ (x,t) is then given by the average of this overlap over the 
ensemble specified earlier. As all collisions have a time- reversed counterpart, the phase of 
the scattering matrix does not matter here and the overlap we are interested in equals the 
two-point correlation function of the classical observable 

g z (x,t) = ^m k 5(x - x k (t)) ; (3.15) 

k 

where we are labelling particles consecutively from left to right with an index k, Xk(t) denotes 
the position of the k particle at time t, and m k is the z component of the spin of the k th 
particle. This correlation function is calculated using the ensemble of initial conditions 
outlined above. The dynamics governing the time evolution of the x k is just that of particles 
moving ballistically except for elastic collisions in which they retain their spin values. 
Thus we can write 

C UjZZ (x,t) =J2( m kS(x - x k (t))mi5(xi(0))) - (g z ) 2 , (3.16) 

k,l 

here the angular brackets refer to averaging over the ensemble of spin labels m^, initial 
velocities v k (0), and initial positions x k (0) specified earlier. Now as the spin-projections 
rrik are not correlated with the initial positions or velocities, the averages factorize. The 
correlators of the m k are easily evaluated as: 
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(m k mi) = A x + A 2 5 k i , 



(3.17) 



where A\ = (f\ — f-i) 2 and A 2 = f\ + /_i — (}\ — f-i) 2 are simple, dimensionless, known 
functions of H/T only. Using ( |3.17| ) we have 



C u , zz (x,t) = A 1 [(p(x,t)p(0,0))- p 2 



(3.18) 



where p(x,t) = Y^k^i x ~ x k(t)) is the spacetime dependent total density, all averages are 
now with respect to initial positions and velocities, and p = (p(x,t)) = J2mPm- The two- 
point correlators of p(x, t) are also easy to evaluate: if the spin labels are neglected, the 
collisions have no effect and correlators of the total denstiy can be obtained by considering 



an ideal gas of point particles. The second correlator in ( 3.18 ), multiplying A 2 , is more 
difficult: it involves the self two-point correlation of a given particle k, which follows a 
complicated trajectory. Fortunately, precisely this correlator was considered three decades 
ago by Jepsen |S| and a little later by others p0| ; they showed that, at sufficiently long 
times, this correlator has a Brownian motion form. In Appendix C, we give a self- contained 
summary of Jepsen's calculation. Here we just write down the final results [22"] for the 
correlation function: 



C u , zz (x, t) — p 2 



L T . Li 



where p 2 F\ is the connected density correlator of a classical ideal gas in d = 1, 

F 1 {x,t) = e~ i2 ^/tV^ , 
and p 2 F 2 is the correlator of a given labeled particle |j49,[50l , 



(3.19) 



(3.20) 



F 2 (x,t) 



(2Gi(«)Gi(--u)+Fi(£,t) 



x/ (2t^G 2 {u)G 2 {-u) 

G 2 (u)G 2 (-u) + G 2 (-u)G 2 (u) 



^G 2 (u)G 2 (-u) 



xh[2tJG 2 (u)G 2 (-u] 



-(G 2 (u)+G 2 (-u))t 



(3.21) 



with u = x/t, Gi(u) = erfc(w)/2, and G 2 (u) = e~ u - uG^u). For \t\ < \x\ < 1, 

the function F 2 has the ballistic form F 2 (x, i) ~ Fi(x, i), while for \t\ 3> 1, \x\ it crosses over 
to the diffusive form 



F 2 (x,t) 



(47rt 2 )V4 



for large t 



(3.22) 
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In the original dimensionful units, ( |3.22j ) implies a spin diffusion constant, D s , given exactly 
by 

c 2 e A/T 

Ds = A(l + 2cosh(Jf/T)) • (3 ' 23) 

This result has been obtained by the solution of a classical model which possesses an infinite 
number of local conservation laws: in Appendix A, we explicitly show how the existence 
of these local conservation laws is not incompatible with diffusive spin dynamics. It must 
be noted that the result ( |3.23| ) does not imply that we have rigorously established that the 
ultimate long-time corelations of the quantum model are also diffusive: the reasons for this 
and related comments were made earlier in Section [I] below (|1.5| ) . 

Let us now summarize the calculation of the correlator of the transverse components of 

the magnetization density. The semiclassical prescription for evaluating C U) ^(x,t) is again 

quite straightforward: We begin with an initial state chosen from the same ensemble as 
before. L + (0) acting on the initial state gives zero unless there is a particle at x = with 
spin label m — 0, — 1, in which case it raises the m value of that particle by 1 and multiplies 
the state by a factor of \J~2 (coming from the usual properties of raising operators for the 
spin-one representation of the angular momentum algebra). The resulting state is then 
propogated forward in time with all the particles moving along their classical trajectories as 
before. At time t, the operator L_(x) acting on this state gives zero unless there is a particle 
at x with spin label m = 0, 1, in which case it lowers the spin value of that particle by 1 and 
again multiplies the state by a factor of \/2. This state is then propogated backward in time 
and its overlap with the initial state calculated. C u - + (x, t) is given by this overlap averaged 
over the ensemble of initial conditions. Here, as before, the phase factor of —1 coming from 
each collision does not matter as each collision has a time reversed counterpart. Also, it is 
easy to see that in this case the overlap with the initial state is zero unless L_(x) lowers 
the spin of precisely the particle whose spin was raised by L + (0). Lastly, we see that there 
is an overall factor of e +lHt coming from the unitary time evolution as the total spin of the 
state during its evolution forward in time is greater than the total spin during its evolution 
backward in time by precisely one. Similar considerations apply to C Uj+ _. Putting all of 
this together we see that 

C Ui ^(x, t) = 2p 2 e± lHt A T F 2 (3.24) 

where A T = f + f Tl . 

Now, we may express the NMR relaxation rate in terms of the correlation functions of 
the conserved magnetization density as 

7fr = J2 A *f$ A ^ S^{u N ), (3.25) 

J l a=x,y p^= x ,y,z 

where the local dynamic structure factor Sp^uisf) is defined as 

Sfriutf) = I dte^ Nt C u ^(0,t) ; (3.26) 
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note that we have neglected the q dependence of the hyperfine couplings and ignored the 
contribution of the antiferromagnetic spin fluctuations to the integral over q in ( p. 13 ). At 



this point we have to address an important subtlety that arises in calculating the local dy- 
namic structure factor from the autocorrelation function. We are treating the spin dynamics 
semiclassically to arrive at our expressions for the correlation functions. This gives rise to a 
characteristic 1/t divergence at short times in the corresponding autocorrelation functions. 
This is basically a signature of classical ballistic spin transport; at these short time scales 
collisions play no role. As a result, the integral as written is logarithmically divergent at 
short times. Our semiclassical expressions for the correlation functions do not make sense 
for very short times. This is natural as our whole approach has been geared towards calcu- 
lating these correlations at time scales much larger than 1/T and length scales much larger 
than the thermal de broglie wavelength; our method fails when both these conditions are 



simultaneously violated | 23[| . The semiclassical expressions for C u (0,t) are thus only valid 
for t > e t where e t is a short time cutoff ~ 1/T. Introducing this short time cutoff will give a 
well- defined result for S a p{uj^) at the price of introducing an arbitrary scale e t ~ 1/T; this 
does not seem very promising as our results for S zz {ujn) (S± t (ljn)) will depend sensitively 
upon e t except for very small fields such that we are in the collision dominated diffusive 
regime: JnH C (H <C V-^t)- Note that the range of H for which the results are 

insensitive to the cutoff differs for the transverse components of the local dynamic structure 
factor because of the overall oscillatory factor of e * in the corresponding autocorrelation 
functions (this factor always dominates as 7jv <C 1). However, we can still use our approach 
to compute the S a p(uJis[). The point is that, at very short times, the collisions between the 
thermally excited particles do not matter, and the spin dynamics is ballistic. This means 
that S a p{uj^), for high frequencies (such that un is much larger than the mean collision 
rate ~ ^/L t ), may be calculated exactly by doing a full quantum calculation for a gas of non- 
interacting spin-one particles f3l| . Now, we can expand our semiclassical result (obtained 
by using a cutoff e t ) for uj^ ^> l/L t and match the leading term in this large H expansion 
with the small H limit (for S zz (S± T ) this would be the regime 7 at if <C T (H <C T)) of the 
quantum calculation of Ref |3l| . This, then, will uniquely fix e t and give us results for the 
S a /3(oo^f) that will work reasonably well even for H ~ T (though strictly speaking they are 
valid only in the range JnH <C T (H <T) for S zz (S^)). 

To see explicitly how this procedure works, consider S zz (ujn). It is quite easy to see that 
the ujn \/L t limit of the semiclassical S zz {ujn) is: 

A ' (e H l T + e- H l T )\J^—\ , 



7TC 2 V Ct^N . 



where 7 ~ 0.577216 is Euler's constant. The ujn -C T limit of the full quantum calculation 



reads [31 



^(e H/T + e~ H/T In . 

7TC 2 \ U N J 

Thus we can set e t = 1/4T to match the two logarithms. It is easy to check that the 
same choice works for the transverse correlators. It is now quite straightforward to do the t 



integrals and obtain the following results G2| for the local dynamic structure factor 
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S„(u N )= [A x {ln(TL t ) + ^(V^Mh)} 

+A 2 {\n{TL t ) + <$> 2 {^\u N \L t )]\ , 

S t± (uj n )= ^^J^ {HTL t ) + $ 2 (v^K ± H\L t )} . (3.27) 



The ln(TXt) terms logarithmically violate the purely classical, reduced scaling forms ||36|| . 
and were fixed using the matching procedure just outlined. The scaling functions $i )2 (f2) 



were determined in Ref 22] to be 



K } K 1 4 v / fi(v / 4TTF+2) 1 /2 

(l + fi 2 /^ 2 (fi))V2( 1 + ^(Q)) 

- ln ™ ; 



2^ • ( 3 - 28 ) 



where 7 = 0.57721... is Euler's constant, and = (fi^l + fi 2 /4 - fi 2 / 2 ) 1/2 - Note 

that the above expression for $ 2 (fi) clearly shows the expected crossover from the large 
frequency ballistic behavior <3? 2 (f2 ~~ > 00) = ln(l/0), to the small frequency diffusive form 
$ 2 (Q^0) =tt/(2VT2). 

Let us now use all of this to make contact with the experimental results of Ref [^UJ. For 
this particular experimental setup, the expression for 1/Ti simplifies and to a very good 
approximation we can write j30] 

^ = TiX S xx (u N ) ; (3.29) 



here the relevant hyperfine coupling constant is known |3(J to have the value Ti w (7.5 x 
10 5 ) sec^ 1 (note that we have used units such that h = ks = 1 in our computation 
of the correlation functions and thus time is being measured in inverse Kelvins). To begin 
with, we straightforwardly attempt to fit the field dependent 1/Ti with our results. We use 
the values A = 320 K and c = 3.32A (we are working in units where the lattice constant 
a is set to one) extracted from the susceptibility data |52|. In actual fact, we introduce an 



additional, field- independent background rate Rb that we add on to our theoretical result 
for \/T\. This serves as our fitting parameter; we choose it at each temperature to achieve 
the best agreement with the results of Ref [Q. We show the resulting fits for T = 320, 220, 
and 160 K in Fig [L6|. We see that the theoretical curves account for the field dependence 
of 1/Ti extremely well in this temperature range (of course the agreement for T = 320 K 
should not be taken too seriously as our theory is valid only for temperatures smaller than 
the gap). In particular, the data seems to clearly exhibit the theoretically predicted \j\J~H. 
divergence at low fields which is a characteristic of diffusive spin dynamics. In Fig [17] and 
Fig ITS], we compare the theoretical predictions with the experimental data at T = 120, 



100, 90, 80, 70, and 60 K. At these lower temperatures this divergence seems to get cut 
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off below some threshold field and the quality of the fit deteriorates rapidly. This indicates 
the presence of some spin-dissipation mechanism which becomes significant at these lower 
temperatures and rounds off the diffusive divergence in the local dynamic structure 

factor. Both inter-chain coupling and single-ion anisotropy of the intra-chain coupling are 
expected to contribute to the spin dissipation rate. However, we do not have any real theory 
that can work out the effects of these terms in the Hamiltonian on the field and temperature 
dependence of 1/T\. 

We can only attempt to phenomenologically introduce some spin dissipation in our the- 
oretical results for the spin correlators. Following POJ, we do this by simply introducing 



an exponential cutoff to the long-time tail of the autocorrelation function; thus we write 
C' uxx (0,t) = e~ li: C UtXX (0,t). It is straightforward, though somewhat tedious to work out the 
corresponding local dynamic structure factor by doing the Fourier transform and we will 
spare the reader the details. This now gives us a phenomenological result for \jT\ with an 
additional tunable parameter 7. We choose this spin-dissipation rate at each temperature to 
achieve the best fit with the data. The resulting curves are shown in Fig [19] and Fig ^0| for a 
few representative temperature values. We see that it is possible to fit the data moderately 
well; discrepencies are however clearly visible and it is not clear how much significance to 
attach to the sharp increase in 7 as the temperature is lowered. The quality of our fit seems 
at first sight to be much worse than the corresponding fit to a purely classical diffusive form 
employed in Ref [30]. However, it is important to note that the phenomenological model 
of Ref |3(| used the diffusion constant as an additional fitting parameter; we do not have 
any such freedom. Moreover, both the diffusion constant and the constant background rate 
extracted from the fit in Ref []3(J take on unphysical values below about 100 K [30]. This 
is because, at these lower temperatures, we are in the ballistic regime of spin-transport for 
a significant portion of the H axis and the contribution from the 'free-boson logarithms' 
cannot be neglected. As the crossover to the ballistic regime is already incorporated in our 
form, the present results for the background rate do not suffer from any such obvious prob- 
lems (the diffusion constant of course is just given by (|3.23|) in our approach). In Fig |21] and 
Fig [22| we plot the corresponding values of the spin-dissipation rate 7 and the background 
rate -R& as a function of temperature. The spin dissipation rate is seen to increase rapidly as 
the temperature is decreased. On the other hand, we see that the temperature dependence 
of Rb may be fit approximately by an activated form with activation gap close to 3 A/2. 



IV. HIGH T REGION (T > A) OF THE CONTINUUM a MODEL 

We consider here the possibility that it may be possible to find gapped spin chains 
which satisfy A < J, where J is a typical exchange constant. In this case, it becomes 
possible to access a higher temperature regime where a continuum field theory description 
is possible in the regime A ^ T ^ J. In particular, we expect that the continuum cr-model 
to apply in such a regime [|35|] . It is our hope that such a universal high T regime can be 



experimentally accessed in S — 2 spin chains [p^] . Moreover, the study of such a high T 
regime is of importance as matching its results with the T -C A theory can, in principle, 
help us estimate the values of T to which the low T results can be applied. 



An important property of this regime [35] is that equal-time two-point correlator of n 
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C(x, 0) (Eqn ( |3.3| )) decays at large x with a correlation length 

e~|ln(T/A). (4.1) 

We will shortly determine the exact values of the prefactor and the argument of the logarithm 
in ( |4.1|) . At distances of order or shorter than this correlation length we may crudely expect 
that a weak-coupling, spin-wave picture will hold, and excitations will have energy of order 
or smaller than c£ _1 , which is logarithmically smaller than the thermal energy T; in other 
words 

^ < 1 (4.2) 

T ln(T/A) K ' 

So the occupation number of these spin-wave modes will then be 

1 T 

> 1 (4.3) 



The last occupation number is precisely that appearing in a classical description of thermally 
excited spin waves, which suggests that a classical wave description should yield an appro- 
priate picture of the dynamics of this high T region. However, notice that classical thermal 
effects are only logarithmically preferred, and any predictions of a classical dynamical theory 
will only be correct to leading logarithms. 

We begin our analysis by first focussing on the equal-time correlations in this region. We 
shall use a method originally introduced by Luscher The main idea of Luscher is to 

develop an effective action for only the zero Matsubara frequency (u n = 0) components of n a 
after integrating out all the u n ^ modes (the uo n = modes are related to the equal-time 
correlations via the fluctuations-dissipation theorem and the Kramers-Kronig relations [0). 
This is expected to yield the following partition function for a r-independent field n a (x): 

VnJr)S(,r - Uoxp ( - {N ~ ^ [dx ( ) j ( J. I) 



4 J \ dx 

We have now generalized to a field n a with iV components, and will quote many of our results 
for general N; the physical case if of course N = 3. The coupling constant in (|4.4j ) is written 
in a form such that £ is the exact correlation length: this follows from the easily computable 
exact correlations of Z by interpreting it as the quantum mechanics of a single quantum 
rotor. The value of £ can be computed in a perturbation theory in g on the quantum model 



J]): the tu n 7^ modes can be integrated out using a now standard approach |47 
(N-1)CT c 2lAT o . fdk 
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The integral on the right-hand-side is not ultraviolet convergent. We evaluate it using 
the renormalization procedure discussed by Brezin and Zinn- Justin We introduce a 

momentum scale fi at which coupling constants are defined, and generalize fl2.1|) to a model 
in d spatial dimensions. We now define the renormalized dimensionless coupling 
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9r(p) = tfZig, ( 4 -6) 
where e = 1 — d, and the renormalization constant Z\ is determined in dimensional regular- 



ization to be H55 



Z 1 = l 



(N - 2) 
2tt e 



+ . . . 



(4.7) 



We now need to express ( |4.5| ) in terms of g R , and evaluate the integral on the right hand 
side in d = 1 — e dimensions. Let us display a few steps of the latter evaluation: 



r E 



1 



(2vr) 1 - e c 2 P + u 2 n 



d}~ e k 



(2< 



l-e 



r E 



1 



c 2 &; 2 + uo 2 



duo 



1 



27T C 2 F + LU 2 + T 2 



+c 



l-e 



d 2 " e p 



1 



(2tt) 1 - £ cV + T 2 



1 /T 

c V c 



d l ~ e k 



1 



r( e /2) 

(4vr) 1 - £ / 2 



1 k 1 

— coth — — ■ , 

2k 2 k 2 2VWTT 



(4.8) 



We are only interested in the poles in e and the accompanying constants, and to this accuracy 
the first integral on the right hand side can be evaluated directly at e = 0, while the T 
function yields a pole. Now combining ( |4.6|) , (|4.7| ) and (|4.8|) into ( |4.5|) we find that the poles 
in e cancel (as they must), and 



(N - 1)£T 1 (N — 2) 



\n(cfi/TVg) 



(4.9) 



2c g R (fj) 2tt 

where the constant Q is 

Q = 47re" 7 = 7.055507955 ... . (4.10) 

with 7 Euler's constant. Now we use the conventional relationship between fi and the 
renormalization group invariant Aj^ [55,53(1 
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2tt 



-l/(N-2) 



2tt 



cxp 



(N-2)g R (n) t 



(4.11) 



to eliminate the scale \i from ( |4.9| ). As expected, the coupling gn{p) drops out of the resulting 
expression, and we get 



c{N - 2) 
Ttt(N — 1) 



In 



QT 
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(iV - 2) cA^ 



(4.12) 
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Finally, we can express this in terms of the T = gap A by using the relationship between 
A^jg and A obtained using the Bethe ansatz solution of the cr-model [56| 



A 



i/e) 



l/(N-2) 



cA 



(4.13) 



MS 



r(i + i/(Jv-2))' 

The results ( |4.12|,|4.13D lead to the N = 3 result for £ quoted earlier in (|1.8|) . 

Having obtained the classical model ( |4.4| ) for the equal-time correlations, and the precise 
value of the coupling £ in ( 4.12f) , we now turn to an examination of unequal time correlations 
in the high T region T ^> A. We employ an approach related to that used in the study of 
the quantum cx-model in d = 2 in Ref |57j in a low T region; unlike (|4.4|) , the equal time 
correlations in d — 2 were described by a theory that was not ultraviolet finite, and this 
will lead to significant differences in the analysis and physical properties here. To obtain 
classical equations of motions we clearly need to extend the classical Hamiltonian in ( |4.4|) 
by including a kinetic energy term, expressed in terms of a canonical conjugate momentum 
to n a . The obvious approach is to take the quantum equations of motion, and to simply 
treat the variables as c-number classical degrees of freedom. In particular, we treat the 
rotor-angular momentum L as a classical variable, and augment the classical Hamiltonian 
by the kinetic energy of rotational motion. The moment of inertia of the rotor is related 
to the response of the system to a magnetic field H, and we therefore need to study the 
behavior of Xu i n the T ^> A regime. 

We will determine Xu by strategy similar to that employed above in the computation of 
£: first integrate out the non-zero frequency modes, and then perform the average over the 
zero frequency fluctuations. We choose an H which rotates n a in the 1-2 plane, and define 



n, 



N-l 

t (x, r) = \Jl - tT 2 (x, r)n a (x) + 2J n a (x, r)e aa (x) 



0=1 



(4.14) 



where n a (x), e aa (x) are a set of A" mutually orthogonal vectors in spacetime, and n a (x, r) 
represent the finite frequency degrees of freedom which must be integrated out. We expand 
the partition function to quadratic order in H, drop all terms proportional to the spatial 
gradients of n a (x) or e aa (x) (these can be shown to be logarithmically subdominant to the 
terms kept), and find that the H dependent terms in the free energy density are 



H 2 



2cg 



[n\ + n 2 2 ) I 1 - X^a) ) + ( e ai e w + e a2 e b2 ) (vr a 7r b ) 

V a ) ab 



— y2{e a ie b 2 ~ e a2 e bl )(e cl e d2 - e c2 e dl ) f dxdr (7r a <9 T 7r 6 (0, 0); 7r c <9 r 7r d (:r, r)) 



(4.15) 



C 9 abed 

Evaluating the expectation values of the it fields, and using orthonormality of the vectors 



n a , Zaa, the expression (|4.15|) simplifies to 
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(4.16) 
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Finally to obtain the susceptibility Xu, we have to evaluate the expectation value of the zero 
frequency field n a under the partition function (|4.4j ). This simply yields (nf) = (n 2 .) = 1/N. 
The first frequency summation is precisely the same as that evaluated earlier for £ in ( (4.5| ), 
while the second is explicitly finite in d = 1 and can directly evaluated; in this manner we 
obtain our final result for Xu- 
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(N - 1)T£ (N - 2) 



2c 2 



2vrc 



(N-2) 
Nnc 



In 



' QT * 

A~MS e , 



(4.17) 



We have omitted the form of the subleading logarithms, which are the same as those in 
( [4.12j ). This result was quoted earlier in ( |1.9| ). 

We have now assembled all the information necessary to describe the effective classical 
dynamics in the region T ^> A. The classical partition function is given by the follow- 
ing phase-space functional integral, which generalizes ( |4.4j ) (and we will now specialize the 
remainder of the discussion to the special case iV = 3): 



Z = J Vn a (x)VL a (x)5(n^ — l)5(L a n a ) exp 



He 
' T 



H c = - I dx 




-Lit 



(4.18) 



where L a is the classical angular momentum density, and L a , n a are classical commuting 
variables. The second term in 7i c was absent in ( ^4.4|) , and represents the kinetic energy of 
the classical rotors: integrating out L a we obtain ( |4.4| ). By evaluating linear response to a 
field under which 



H c — J dxH a L a . 



(4.19) 



we find 



Xu 



N 



XuJ 



(4.20) 



with N = 3 (we have given, without proof, the expression for general N); the factor of 2/3 
comes from the contraint L a ■ n a = 0. Using (|4.17| ), we then have the value of Xu±- 

We can finally specify the manner in which time-dependent correlations have to be com- 
puted in this effective classical model. The classical equations of motion are the Hamilton- 
Jacobi equations of the Hamiltonian 7i c , with Poisson brackets which are the continuum 
classical limit of the quantum commutation relations : 



{L a (x),L f3 (x')} PB 
{L a (x) } np(x')} PB 
{n a (x),nf3(x')} PB 



e^L^^S^x — x) 
e aj g 7 n 7 (x)5(x — x') 
0. 



(4.21) 



The equations of motion are 



40 



dn a 

dL a 
dt 



^XuA 



(TO 



(4.22) 



The classical correlation functions are obtained by averaging these deterministic equations 
over an ensemble of initial conditions specified by fl4.18|) . Note also that simple dimensional 
analysis of the differential equations ( f4.22| ) shows that disturbances travel with a character- 
istic velocity c(T) given by 



c(T) = {TiiT)l Xul _{T)fl\ 



(4.23) 



Notice from ( |4.12j ) and ( [4.17| ) that to leading logarithms c(T) pa c, but the second term in 
the first equation of ( |4.17| ) already shows that exact equality does not hold. 

We complete the relationship of the quantum to the classical model, by noting that there 
is also an additional wave-function renormalization of the n a field |47]j55| which appears when 
the non-zero frequency modes are integrated out. Our final result for the correlator C in 
73|) then takes the form 



C a p{x,t) = AG 



In 



(iv-i) 

(N-2) 



(4.24) 



The subscript c represents the classical average specified by ( |4.18|) and ( |4.22| ). The constant 
A is the T = quasi-particle residue which appeared in ( |3.5|) . The constant Q is an 
unknown universal number which cannot be obtained by the present methods. It could, in 
principle, be obtained from the Bethe-ansatz solution. There is no similar renormalization 
of the correlator of the magnetization density , C u in (|3.14j ), which is precisely equal to the 
two-point correlator of L a under ( 4.18 ) and (|4.22j ). 

It is now possible to perform a simple rescaling and to show that the classical dynamics 
problem above is free of any dimensionless couplings, and is a unique, parameter-free theory. 
This will allow us to completely specify the T dependence of observables upto unknown 
numerical constants. Let us perform the following rescalings on ( |4.18j ) and ( |4.22| ) 



X = Xt; 




(4.25) 



Then the partition function ( |4.18| ) is transformed to 



Z : 



T>n a (x)VL a (x)5(n 2 a — l)5(L a ■ n a ) exp {—Hi 

2 



dx 



dn a 
dx 



(4.26) 
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while the equations of motion become 



dn a 

dLa 
dt 



d 2 n^ 



dx 2 



(4.27) 



Notice that coupling constants and parameters have been scaled away, and ( 4.26| , 4.27 ) con- 
stitute a unique problem that must be solved exactly. The T and A dependencies of all 
quantities arise only through the rescalings defined in (|4.25|) and the results ( |4.12j ) and 
( [4.17||r2"UD for £ and Xu± given earlier. Complete description of the correlators now requires 
exact solution of ( |4.26|J4.2^) . The equal-time correlations are of-course known from ( f4.26| ): 



L a (lS,0)Lp(0,0)) z = -6 a pS(jE) 
(n«(3f, 0)^(0, 0))_ = l -5 aP e-W 



(4.28) 



Even though the equations of motion constitute an integrable system with an infinite number 
of non-local conservation laws p)8"| , |5TSfl , it is not known how to analytically compute correla- 
tions averaged over the initial conditions of a thermal ensemble, or whether the correlator 
(L a (x,t)Lj3(0, 0))c has a diffusive form at long times and distances. If diffusion did exist in 
the continuum equations ( |4.28j ), the present analysis does allows us to completely specify 
the T dependence of the diffusion constant; by a simple dimensional analysis of ft4.25|) , we 
get 



D„ = B 



TV 2 [£(T)] 



3/2 



(4.29) 



where B is an unknown universal number, and the T dependencies of £ and Xu± are in ( f4.12| ) 
and (|4T7|J4T20| ). 

In this context, it is interesting to note that recent measurements [BD] of the field depen- 
dence of 1/Ti in the compound (VO)2P20y at temperatures T ^> A seem to provide clear 
evidence for spin diffusion. However, the bulk of the data is at temperatures comparable 
to the microscopic exchange constants of the system and it is not clear if the foregoing 
description based on the universal high temperature properties of the continuum field the- 
ory is applicable in the temperature regime studied experimentally. It is interesting that 
the experimental results appear to suggest that D s ~ c£, which is consistent with ( |4.29| ) if 
Xu ~ T£/c 2 (as is the case with our results (|1.8j) and (|1.9|) to leading logarithms). 



V. CONCLUSIONS 

The main results of the paper are already summarized in Section |, and here we will note 
some unresolved issues and directions for future work. 

All experimental realizations of gapped antiferromagnets have additional complications 
which have not been included in the model systems studied here. Most important among 
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these are the spin anisotropies away from perfect Heisenberg symmetry and the inter-chain 
couplings which make the system only quasi-one-dimensional. 

Consider first the consequence of anisotropy. The three-fold degenerate quasiparticle 
spectrum will now be lifted, and three resulting particles will have have different energy gaps 
and masses. Further, these parameters will depend in a complicated way upon the external 
field. Nevertheless, we expect that the simple structure of the iS-matrix in Q1.2Q will be 
retained, as it only depends upon simple dimensional properties of slowly moving particles 
with a quadratic dispersion. Correlations of the particle density can probably be computed 
along the semiclassical lines of this paper: one has to deal with a classical gas of particles 
of different masses and average densities. The latter problem is considerably more complex 
than the equal mass case, and there is probably no alternative to numerical simulations. 
Correlations of the spin operators appear more problematical- these will invariable change 
the labels of the particles when they act, and therefore lead to differences in the labels in 
the forward and backward trajectories. Combined with the complication that the masses 
of the different particles are different, and so their trajectories will have different velocities, 
we are faced with what appears to be a very complex problem with quantum and classical 
effects intertwined. 

Inter-chain couplings will eventually require us to consider dynamics in two or three 
dimensions. If temperatures are low enough that the inter-chain motion is coherent, then 
we have to consider the 5-matrix for scattering in higher dimensions. In this case the 
low- momentum behavior is quite different: in fact the T-matrix vanishes at low momenta 
for d > 2. We would then expect all scattering to be dominated by elastic scattering of 
impurities which would control the behavior of the spin diffusion constant and the quasi- 
particle broadening. On the other hand, systems with only incoherent hopping between 
chains will probably be dominated by the inelastic scattering along the one-dimensional 
chains, and display behavior qualitatively similar to that discussed in this paper. 
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APPENDIX A: LOCAL CONSERVATION LAWS AND SPIN DIFFUSION 



The computation of the spin diffusivity in Section |III B| was carried out using the exact 
solution a simple classical model of point particles in one dimension. This model is exactly 
solvable |^9| and possesses an infinite number of local conservation laws, as we will show 
explicitly below. The existence of spin diffusion then appears to run counter to the conven- 
tional wisdom that the time evolution of a integrable system is not 'chaotic' enough to be 
compatible with diffusion. In particular, one might expect that any non-zero spin current 
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produced in the system will not ultimately decay to zero because the numerous conserva- 
tion laws prevent it. In this appendix we will show that this expectation does not hold 
for the particular model being studied, and that an important 'particle-hole'-like symmetry 
allows complete decay of any spin current at H = 0. In a finite magnetic field (H ^ 0), 
the particle- hole symmetry is absent, and then the spin current does not decay completely: 
this is consistent with the presence of a purely ballistic component, Fi, in ( |3 . 1 9|) which con- 
tributes only for H ^ (A\ = at H = 0), and the arguments of Zotos et al jpU| . A closely 
related particle-hole symmetry also played an important role in the appearance of a finite 
conductivity in our recent quantum transport analysis in two dimensions pT| . 

The classical model of Sections |111 A| and [U1B| consisted of particles k — 1 . . . N with spins 
nik chosen randomly (at H = 0) from 1, 0, —1. At time t = the particles had uncorrelated 
random positions Xfc(0), and subsequently they occupied 'trajectories' Xk(t) = £fc(0) + v k t 
where v k are uncorrelated random velocities chosen from a Boltzmann distribution. The 
position Xk{t) of particle k was however a rather complicated function of time, and was 
chosen from the set of trajectories, {X/~(t)}, such that for all t, Xk{t) < xi(t) for every k < I. 

It is useful at this point to note two discrete symmetries of the above classical statistical 
problem at H = 0. The first is the time-reversal symmetry, T, under which both spins and 
velocities change sign: 



v k -> -v k 



m k 



-m k . 



(Al) 



The second is the 'particle-hole' symmetry P, under which only the spins reverse direction: 



P : v k -> v k 



-m k . 



(A2) 



These symmetries will be crucial in our discussion below. 

Let us now explicitly identify the local conserved quantities of this classical dynamics. 
All of the velocities v k are clearly constants of the motion. However, we would like to work 
with locally conserved quantities which can be written as the spatial integrals over local 
observables, and which are invariant under permutation of the particle labels; so we define 



dx 



N 

£' 

k=l 



(A3) 



with n = 1 . . . N (notice dxk{t)/dt ^ v k = dX k (t)/dt, but the result holds after summation 
over k because the set {x k (t)} differs from the set {X k (t)} only by a renumbering). All the 
V n are constants of the motion. Similarly, with spins mk we can define 



N 

E 

k=l 



HI, 



(A4) 



with p = 1, 2, as additional locally conserved quantities (Mi is the spatial integral of g z (x, t) 
in ( |3.15|) , and a similar result holds of M 2 ). We can now easily work out the signature of 
the V n and M p under the discrete symmetries noted earlier, and tabulate the results: 



44 





P 


T 


V n , n odd 


1 


-1 


V n , n even 


1 


1 


Mi 


-1 


-1 


M 2 


1 


1 



(A5) 



The central quantity in spin transport is the total spin current J(t), which is not a 
constant of the motion. It is also given by an integral over a local quantity as 



J(t) 



dx 



J2 mk ~^~ ^ x ~ Xk ^ 



N 



mjV k Aj k (t), 

j,k=i 



(A6) 



where Ajk is defined to be equal to 1 if particle j is on trajectory k at time t and otherwise; 
we will analytically study the function Ajk in Appendix C, but here we will be satisfied by 
a numerical simulation. Again, as in (|A5|) it is useful to note the signature of J under the 
discrete symmetries: 



(A7) 



As will become clear shortly, one of the central points of this appendix is that the signatures 
in ( |A7| ) differ from all of those of the conserved quantities in ( |A5| ). The current J(t) is the 
sum of iV random numbers of each sign, and so is expected to be of order y/N for a typical 
initial condition chosen from the ensemble defined above. We show the deterministic time 



evolution of J(t) for one such initial condition for a system of 400 particles in Fig 23: notice 
that it is rather noisy-looking and repeatedly changes its sign. Also, among the constants 
of the motion above, we expect V n with n odd and M\ to be of order \/N (provided n is 
not too large), and V n with n even and M 2 to be of order N for a typical initial condition; 
notice that it is only the conserved quantities of order yN that can distinguish left movers 
from right movers, or spin up from down. 

Let us now create a macroscopic spin current (of order N) in this system. We do this by 
hitting the system with a magnetic field gradient impulse at a time t = to, and subsequently 
setting the field to zero. As a result of the impulse, the velocities of the particles with spin 
up are assumed to increase by v , while those of spin down are assumed to decrease by v . 
Formally, this can be written as 



v k — ¥ + rriiVo where / is unique solution of A lk (t ) 
Immediately after the impulse, J(t) will have a macroscopic value 

J(t+) = ^Nv + 0{VN) 



1. 



(A8) 



(A9) 



The subsequent deterministic time evolution of J(t) is also shown in Fig it decays in a 
few collision times to a value of order yN and then appears to chaotically oscillate in time! 
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The basic point is now easy to see. Because is as likely to be +1 or —1, the impulse 
on any given particle is equally likely to be +v or — v . Hence the V n , with n odd, remain 
of order even after the impulse. This is simply a manifestation of the fact that the 
signatures of J under P and T are different from those of the V n . A non-zero J is therefore 
not correlated with an induced value of a conserved quantity which could prevent the decay 
of J to non-macroscopic values. 



APPENDIX B: NUMERICAL COMPUTATION OF THE FOURIER 
TRANSFORM OF THE CORRELATION FUNCTION C 

In this appendix we outline the numerical method employed in calculating S(q,u) start- 
ing from the numerically determined semiclassical C (x, t) and the procedure used to directly 
determine the scaling function (see Eqn ( |3.12|) ) 

As the numerical determination of R(x, t) is the most time consuming part of the entire 
procedure, we calculated R only at a predetermined grid of points in the x — t plane. We 
chose t values from to 7.0 at intervals of 0.2. For each such value of t, we chose about 
20 points so as to sample R as well as possible in the region in which R > 5 x 10~ 3 ; this 
choice was made to reflect the fact that our absolute error in R was estimated to be about 
5 x 10~ 4 . This then defined our grid. At each t, we fit R as a function of x to the form 



, / f,x a i a 2 + CL3X + a A x 2 + f 
\og(R) = — — — 

02 + a 5 x + x z 



where / = 4/3 and a\ = — log(i?(0, t)). The rationale behind our choice of the value of / is as 
follows: When x ^> t, the complicated correlations between the spin labels of a given classical 
trajectory at different times do not matter and R is well approximated by our 'mean-field' 
theory (see Appendix D). The mean-field theory in this limit gives log(i?) ~ — 4x/3 and 
this is what determines our choice of /. The error in the fit was estimated to be roughly the 
same as the error in the original computation of R; thus we did not lose anything by doing 
the fit. Having tabulated the fitting parameters for each value of t on the grid, we evaluated 
the spatial Fourier transform numerically. The resulting function of t is expected to be 
smooth as long as 86j = {u — e(k)) L t is not too large. More precisely, we do not expect any 
oscillations on the scale of our grid spacing in t as long as 0.25d) -C 2-n. As we are interested 
only in 5 Co ~ 1, we can safely interpolate the resulting function in t. In practice we use 
a cubic-spline to do the interpolation. Lastly, we do the t integral numerically to obtain 
S(q,u). The accuracy of both numerical integrations is quite high and so we expect that 
the dominant error in our calculation comes from the interpolation; this is conservatively 
estimated to be a few percent at the most for the largest values of 50. 

Let us now briefly indicate the procedure used in obtaining the Fourier transform of 
-R(0, t) needed for the calculation of the scaling function <&(z). The available data for R(0, t) 
is fit extremely well by the following form: 

. , ~ . at + at 2 + bt 3 
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where the choice a = ^jZ^pn is again motivated by 'mean-field' considerations. It is now 
a simple matter to do the Fourier integral to a very high accuracy using this fit and we 
estimate the errors involved to be less than 0.5% at the most. 



APPENDIX C: CALCULATION OF TAGGED PARTICLE CORRELATIONS IN 

THE CLASSICAL MODEL 

In this appendix, we shall attempt to give a self-contained account of the method devised 



by Jepsen [IS] for the calculation of the tagged particle correlations in the classical model 



introduced in Ref |49 |. We will try to adhere to the notation and conventions of [[49| as far 
as possible. 

The model is defined as follows: We begin with N particles of mass m distributed 
uniformly along a one-dimensional segment of length L with periodic boundary conditions 
(we will eventually take the thermodynamic limit L —>■ oo with N/L fixed to be equal 
to the density p). At time t = each particle is assigned a velocity from the classical 
thermal ensemble defined by the usual Maxwell-Boltzmann distribution function g(v) = 
(m/2nT) 1 ^ 2 e~ mv l 2T . The subsequent evolution of the system is purely deterministic; the 
particles travel without any change in their velocities until they collide with another particle. 
Every collision is elastic and the particles merely exchange their velocities as a result of the 
collision. 

To begin our analysis, let us label the particles from left to right with an index i running 
from to N — 1. Thus the particles are initially at positions 2,(0) such that Xi(0) < Xj(0) 
for i < j. Actually, it is convenient to identify i + N with i because of the periodic boundary 
conditions employed which identify the ends x = and x = L of the interval. Note that this 
labelling of the particles is left invariant by the dynamics. We also label trajectories (which 
follow the straight line defined by Xi(t) = Xi(0) + Vit on the space-time diagram representing 
the evolution of the system) with an index i, again with the convention that Xi(0) < Xj(0) 
for i < j (here v i is the initial velocity of the i th particle). Let Xi(t) denote the position of the 
i th particle at time t. We wish to calculate the correlator B(x,t) = (5(x — Xk(t))5(xk(0))) 
where summation over the repeated index k is implied and the angular brackets refer to 
averaging over the ensemble of initial conditions specified earlier. 

Let us now consider the quantity Ajk(t), introduced in Appendix A, which is defined 
to be equal to 1 if particle j is on trajectory k at time t and otherwise. Another useful 
quantity is the number nu of (signed) crossings suffered by the k th trajectory upto time t. 
Every time this trajectory is hit from the left, decreases by 1 and every time it is hit 
from the right ni~ increases by 1. Clearly, A^it) = 1 for j = k + nk(t) and zero otherwise. 
We may probe the dynamics a bit more by defining another quantity r n (h, k, t) which equals 
1 if trajectory h has crossed trajectory k precisely n times upto time t and zero otherwise. 
Here too, we are talking of signed crossings; if trajectory h crosses from the left this is a 
negative crossing and if it crosses from the right it is a positive crossing. Clearly r n has the 
interpretation of a probability when averaged over any ensemble of initial conditions. Let 
us also define the corresponding 'generating function' as 

oo 

s(u;h,k,t)= r n (h,k,t)e mu . (CI) 
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The reason for introducing r n and s{u) is that Ajk(t), which is clearly a central quantity 
of interest, may be very conveniently expressed in terms of s(u) as 



i N-l 

A A*) = ^ E e ~~ 

iV 1=0 



N-l 

n 

m=0 



'2tt/ 



■; m, k, t 



(C2) 



here we are using the convention that s(u; k, k,t) = 1. This is quite easy to check from the 
defmtions of s(u) and Aj k . Moreover, it is possible to write down a fairly explicit expression 
for s(u; h, k, t). This takes a slightly different form depending on whether h is greater or less 
than k. If h > k, we have 



s(u; h, k,t) = S [u,w kh ] , 

while if h < k, we have 

s(u; h, k, t) = e~ m S [u, w kh ] . 
S [u, Wkh] used above is defined as: 

S[u,w kh ] = e mu , 



(C3) 



(C4) 



(C5) 



where w k h = x k (0) — Xh(0) + (v k — Vh)t, and n is the integer that satisfies (n— 1)L < w k h < nL. 
Using this definition, we can write the following compact expression for Aj k (t) in terms of 
S: 



1 N-l „ , N-l 

= e ^ n ^ 

iV Z=0 h=0 



2tt/ 



(C6) 



With all this machinery in place, it is a relatively straightforward matter to calculate 
the correlation function we need. We begin by explicitly writing out the ensemble averages 
involved: 



B(x,t) = jj 7 J {dx} J [dv] ^ 9(vi)j 5(x (t)-x) 



(C7) 



where we have used the definitions / {dx} = J L dx^-i Jo 1 *' 1 dx^-2 ■ ■ ■ Jo 2 dx\ and / [dv] = 
/^o dvpj-i ■ ■ ■ I^°oo d v o with x k (0) = x k , and it is understood that x (0) is set 

equal to when evaluating the right hand side of (|C7|) . Now we can transform from particle 
positions to trajectories by writing 

5(x (t) - x) = J2 A okS(X k (t) - x) . 

k 

Using this and writing A ok in terms of S [u, w k h] allows us to express our correlation function 

as 



B(x,t) =—j [dx] j [dv] n 9{vi) Ur e n s 

' \ 1=0 / k,l=0 h=0 



2tt/ 



w kh 



S(X k (t)-x). (C8) 
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Here we have also used the fact that the integrand in this representation is explic- 
itly symmetric in the spatial integration variables to change the spatial integration to 
/ [dx] = Jo dx N _i Jq dx N _ 2 ■■■Jo dx 1 . 

It is now convenient to define R [u, x k + v k t — xy\ = /f^ dvhg{vh)S [u, w k h]. Using this 
we can rewrite our expression for the correlation function as 



B(x,t) = -^Yl J dv og(v )5(X (t) -x) (tJ q dx h R(u,X (t) - x h ) 



N-l 



N 



dxi 



dv k g{vk)R{u,Xk{t))6(X k (t) - x) 



x \\ J dx h R(u, X k {t) - x h ) 



N-2 , 



(C9) 



where k ^ h, k,h 0, u = 2irl/N and J2 U = I^Lo 1 - To proceed further we need to work 
out R [u, X k (t)] and (1/L) J L dxhR [u, X k (t) — Xh]- This is quite straightforward to do in the 
limit of large L and we only give the final results below: 



R [u, X k {t)\ = -E c {y) + 1 - -E c (y) e 



Y I (l.r h R[u. X)At) - .(■/,] = 1 + (1 



L Jo 



T[u,X h (t)) 



(CIO) 



where 



2T 

T[u,X k (t)) = \—t 



m 



E c (y) 



y 



ye m + 



dze 



i 1/1 



TC Jy 

rnX k {t) 



2T t 

Now, in the the thermodynamic limit specified earlier we can write 



(Cll) 



l + I (l- e - iu )T[u,X k (t)} 



N-v 



exp(p(l-e- iu )T[u,X k (t)} 



valid for any finite number v. Using this and ( CTD ) in the expression (|CJ]) for the correlation 
function and doing the remaining integrals over positions and velocities gives us 



B(x,t)=p 



27T du 
o 2^ 



1 



+ 



exp 

e~ w2 exp(p(l-e- iu )T[ 
t V 2nT P v v 1 



e~ m )T[u,x 



m 



u,x 



(C12) 



here we have replaced the sum over u by the corresponding integral in the thermodynamic 
limit and used the definitions f\(w) = E c (w)/2, f2(w) = 1 — fi(w) and w = (m/2T) 1 ^ 2 x/t. 
To do the u integral, we note that T may be expressed as Ge tu — A where A and G are 
functions purely of x and t. This allows us to use the standard Bessel function identity, 
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r2n 



2tt Jo 



due 



exp 



L )(Ge iu - A) 



A 



e~ p{A+G) I n (2p^AG) 



(C13) 



to finally arrive at the results quoted in ( |3.21| ) of Section [111B| upon using the appropriate 
values for p and m. 



APPENDIX D: APPROXIMATE ANALYTICAL CALCULATION OF THE 

RELAXATION FUNCTION 

In this appendix, we briefly outline our approximate 'mean-field' theory for the relaxation 
function R(x, t). 

We begin by noting that the classical model defined in Section |III A| has been solved 
exactly in Ref |23| for the special case in which there is only one possible value for the spin 
label m. All of the difficulties we encounter in attempting to generalize this solution to the 
case of interest here stem from the fact that there are complicated correlations between the 
rrik(t) (defined in Section [III A|) at different times. 



Our 'mean-field' approximation consists of simply ignoring these correlation effects 
(hence our choice of terminology to describe our approximation). Having made this un- 
controlled approximation, it is now a fairly straightforward matter to obtain a closed form 



expression for R(x,t) in analogy with the corresponding discussion in J2I|. The actual cal- 
culation proceeds as follows: Let q be the probability that any given solid line in Fig |9| 
intersects the dotted line. If we ignore the correlations between the m^{t) at different times, 
then the probability that this line carries a spin label equal to the spin label of the dotted 
line is 1/3. So given that the line intersects the dotted line, this intersection contributes a 
factor of —1 to R(x, t) with probability 1/3 and a factor of with probability 2/3 (if the line 
does not intersect at all, we of course get a factor of 1). Within our mean field theory, R is 
just a product of such factors, one from each solid line. This gives R(x, t) = (1 — q — q/3) N , 
where N is the total number of thermally excited particles in the system. Now, using 
q = (\x — vt\) / L [^] (where the angular brackets denote averaging over the Maxwell- 
Boltzmann distribution function for v and L is the length of the system) and taking the 
thermodynamic limit, we obtain R(x,t) = exp(— 4p(\x — vt\)/3). We can now do our usual 
rescalings and write down the main result of our mean field theory: 

R(x, t) = exp(-4(|x - i5f|)/3) , (Dl) 

where the angular brackets now denote averaging over the distribution V(v) = -^e~ v2 and 

x and t are defined as in Section [III A| . In particular, note that this implies R(0, t) = 
e -4|t|/3V7r. t urns out to be reasonably accurate for some purposes (see the discussion on 
the approximate form of the scaling function $(z) in Section [111 A] ). 
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FIG. 1. Two particle collision described by the 5-matrix ( |1.2| ). The momenta before and after 
the collision are the same, so the figure also represents the spacetime trajectories of the particles. 
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FIG. 2. Low and high temperature asymptotics for the uniform susceptibility \u of the contin- 
uum 0(3) non-linear cr-model. At T = 0, there is an energy gap A to all excitations, and c is the 
velocity defined by The expression in Eqn 1.4 gives the low temperature asymptotics while 

Eqn |1.9| is used for the high temperature asymptotics. Any lattice antiferromagnet will have a very 
high temperature (T > J where J is a typical microscopic exchange constant) Curie susceptibility 
~ 1/T which is not shown: the high temperature limit of the continuum theory will apply for 
A < T < J. 
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FIG. 3. Positions in u> of the single particle peak (solid line), bound state peak (long-dash line), 
and the bottom of the two particle continuum (short-dash line) in S(Q,uj) plotted as a function of 
q x for the strongly-coupled ladder (a typical value of g = 0.25 is used). 
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FIG. 4. Spectral weight in the single particle peak (dashed line) and the bound state peak 
(solid line) in S(Q, u) for a strongly coupled ladder (a typical value of g = 0.25 is used). Note that 
the two curves actually correspond to different values of the transverve momentum q y chosen to 
maximize the respective spectral weights: the single particle part is shown for q y = ir/d while the 
bound state part is shown for q y = (d is the spacing along the rung of the ladder). 
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Propogator denoted as 
on the left. 

Associated expression 
is denoted as D( k, w). 




Vertex (a) 




Vertex (b) 



FIG. 5. The Feynman rules we need for the calculation described in Section II C. The pro- 
pogator D(k,u) is given as D(k,u) = i/(u — e(k) + iff). The factor corresponding to the 
vertex (a) is ig(fR{pz) /r(pi) /i(P4) Il(P2) + R < — ► -^)/2 • The factor corresponding to (b) is 

i9(fLip^)fR(pi) - fR(jPl)h{p\)){fL{jPz)fR(jP2) ~ fR(P3)h(P2)) • 
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FIG. 6. Feynman diagrams contributing to T^k x,—k x;k y,—k y) to first order in g. All 
external lines carry on-shell frequencies corresponding to the momentum labels shown. 
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FIG. 7. Feynman diagrams contributing to T^{k x,—k y;k y,—k x) (diagram (a)) and 
T/±(k x, —k y; k x, —k y) (diagram (b)) to first order in g. All external lines carry on-shell fre- 
quencies corresponding to the momentum labels shown. 
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FIG. 8. Ladder series giving the leading infrared divergent terms in the expansion for iM% 
(diagram (a)) and 1M2 (diagram (b)). All external lines carry on-shell frequencies corresponding 
to the momentum labels shown. The internal lines also carry frequency labels that are not explicitly 
shown. 
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FIG. 9. A typical set of particle trajectories contributing to C(x,t). Each full line represents 
paths traversed by particles moving both forward and backward in time. The dashed line is 
traversed only going forward in time. Shown on the trajectories are the values of the particle spins 
mfc which are independent of t in the low T limit. Notice that all the trajectories intersecting the 
dashed line have a spin equal to that of the dashed line: only such configurations contribute to the 
C(x,t). 
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FIG. 10. S(q,ui) rescaled by a factor of AcLt / '(tt 2 A) plotted against 56) = Lt(uj — A) with 
q = n/a for A/T = 2, 3, 4, and 5. Note the scaling collapse of the curves corresponding to the 
three lowest temperatures. 
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FIG. 12. The scaling curve of Fig 11 (defined by the data for A/T = 4 and A/T = 5) plotted 
against the independent variable (— fc 2 /2) compared to the scaling curve of Fig |l0| (again defined 
by the data for A/T = 4 and A/T = 5) for the corresponding negative values of 5uj. The two 
coincide within our numerical errors. 
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FIG. 13. S(k + ir/a, uj) rescaled by a factor of AcLt/(ir 2 A) plotted against 5uj = Lt (uj — 2 1//2 A) 
with k = A/c for A/T = 2, 3, 4, and 5. 
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FIG. 16. Field dependence of \/T\ for T > 120 K. The experimental data of [3C] is compared 
with the theoretical predictions offset by a field-independent background rate R^ which is the only 
free parameter of the fit; the fit value of is shown under the theory column. 
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FIG. 17. Field dependence of 1/T\ for a few temperatures T < 120 K. The experimental data 
of [[}(]] is compared with the theoretical predictions offset by a field-independent background rate 
Rb which is the only free parameter of the fit; the fit value of Rb is shown under the theory column. 
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FIG. 18. Field dependence of 1/Ti for the lowest temperatures for which data is available. The 
experimental data of |3(| is compared with the theoretical predictions offset by a field-independent 
background rate Rb which is the only free parameter of the fit; the fit value of Rb is shown under 
the theory column. 
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FIG. 19. Field dependence of 1/T\ fit to the phenomenological form described in the text. The 
experimental data of [30| at T = 120, 100, and 90 K is compared to our phenomenological form 
that incorporates a spin dissipation rate 7 in addition to a field independent background rate Rb- 
The values of Rb, 7 are listed under the theory column. 
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FIG. 20. Field dependence of 1/Ti fit to the phenomenological form described in the text. The 
experimental data of |k| at T = 80, 70, and 60 K is compared to our phenomenological form that 
incorporates a spin dissipation rate 7 in addition to a field independent background rate Rb- The 
values of i?6,7 are listed under the theory column. 
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FIG. 21. Temperature dependence of the spin-dissipation rate 7 determined by fitting our 
phenomenological form for 1/T\ to the experimental data of 
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FIG. 22. Temperature dependence of the background rate Rb determined by fitting our phe- 
nomenological form for 1/T\ to the experimental data of pOfl . We plot ln(Rb) against 1/T to check 
for activated behaviour and indeed find an approximate linear relation, the best fit for the slope 
being 468 K. 
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FIG. 23. Deterministic time evolution of the spin current J(t) (defined in (|A6|)) for two systems 
of 400 particles on a circle with the same initial condiions; the value of J(t) changes in discrete steps 
at each collision between a pair of particles. For one of systems, there is an impulse in velocities 
given by ( |A8| ) at a time to = 2. This produces a macroscopically significant J(t), which however 
decays away in a few collision times. The only remnant of the impulse is a 'heating' of the system, 
reflected in the larger amplitude of the order \/N fluctuations in J(t) for the impacted system. 
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